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FOREWORD 


During  the  past  five  years,  the  specialized  engineering  analysis  techniques 
required  in  solid  propellant  grain  structural  integrity  evaluation  have  received 
increasing  attention.  Following  the  Joint  Services  Physical  Property  Panel 
Meeting  in  1958,  where  the  need  for  research  in  this  area  was  emphasized, 
intensified  effort  resulted  in  development  of  engineering  methods  for  analysis 
of  viscoelastic  propellant  grains.  A  program  of  research  under  the  direction 
of  Dr.  M.  L.  Williams  at  the  California  Institute  of  Technology  resulted  in 
the  publication  of  "Fundamental  Studies  Relating  to  Systems  Analysis  of  Solid 
Propellants"  (GALCIT  SM  61-5)  in  1961.  The  report  rapidly  gained  the  status 
of  a  reference  text  for  the  solid  propellant  industry,  and,  with  the  permission 
of  the  authors,  has  been  excerpted  extensively  in  this  handbook.  Indeed,  the 
present  handbook  is  regarded  as  an  extension  of  the  GALCIT  SM  61-5  report, 
with  primary  attention  directed  at  presentation  of  engineering  methods  in  pro¬ 
pellant  grain  structural  analysis. 


The  physical  behavior  of  solid  propellant,  in  contrast  to  elastic  solids,  is 
significantly  affected  by  small  changes  in  the  temperature  and  characteristic 
time  scale  under  consideration.  For  this  reason,  the  methods  of  linear  visco¬ 
elastic  stress  analysis  are  used  to  quantitatively  assess  the  behavior  of  a 
grain  under  pressure,  temperature  and  environmental  loading.  This  handbook 
does  not  attempt  to  treat  the  complexities  of  viscoelastic  structural  analysis 
theory  in  detail,  but  rather  to  describe  engineering  methods  that  have  been 
found  useful  in  applying  this  theory  in  grain  analysis.  With  exceptions  as 
noted  in  the  text,  the  engineering  methods  described  herein  have  been  verified 
in  motor  structural  models,  static  test  motors,  and  flight  motors  in  programs 
at  Lockheed  Propulsion  Company. 


The  handbook  has  been  compiled  for  engineering  use.  Section  1  describes, 
without  extensive  background,  methods  for  laboratory  characterization  of 
propellants.  Section  2  treats  state-of-the-art  engineering  grain  analysis 
methods,  also  without  extensive  background.  More  specialized  aspects  of 
these  problem  areas  are  discussed  in  Sections  3  and  4  and  pertinent  litera¬ 
ture  source  references  are  given.  The  references  cited  are  mainly  those  to 
which  convenient  access  is  possible;  no  priority  of  publication  is  implied. 
Reference  data,  propellant  specimens  preparation  techniques,  and  auxiliary 
laboratory  methods  are  included  in  the  Appendices.  The  information  presented 
in  this  handbook  is  comprehensive  only  in  that  sufficient  information  has  been 
included  to  permit  complete  engineering  analysis  of  grain  structural  problems 
presently  regarded  as  having  been  reduced  to  practice.  Grain  structural 
problems  associated  with  complex  motor  designs  and  complex  material 
properties  require  highly  specialized  engineering  and  mathematical  treatments 
that  are  considered  beyond  the  present  scope  of  the  handbook. 
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Advantage  haa  been  taken  of  the  fact  that  extensive  information  concerning 
propellant  physical  characterization  will  be  available  about  mid  1963  in  the 
"ICRPG  Solid  Propellant  Mechanical  Behavior  Manual".  Summary  infor¬ 
mation  concerning  certain  test  procedures  has  been  included  herein  with 
reference  to  more  extensive  discussions  in  thelCRPG  Manual,  as  appropriate. 


This  handbook  has  been  compiled  and  edited  under  the  direction  of  J.  W.  Jones 
by  the  Lockheed  Propulsion  Company  Structural  Integrity  Department. 


-v- 


LOCKHIKD  PROPULSION  COMPANY 


Section  1 


ENGINEERING  PROCEDURES  FOR  PROPELLANT 
PHYSICAL  CHARACTERIZATION 

CONTENTS 

Page 

1.  ENGINEERING  PROCEDURES  FOR  PROPELLANT 

PHYSICAL  CHARACTERIZATION .  1-3 

1.1  MASTER  STRESS  RELAXATION  MODULUS .  1-4 

1.1.1  Reduced  Variable  Data  Reduction .  1-9 

1.2  CONSTANT  CROSSHEAD  RATE  ULTIMATE 

PROPERTIES .  1-15 

1.3  CONSTANT  LOAD  ULTIMATE  PROPERTIES .  1-21 

1.4  CONSTANT  STRAIN  ULTIMATE  PROPERTIES .  1-24 

1.5  PROPELLANT-LINER-SUBSTRATE  BONDING 

ULTIMATE  PROPERTIES .  1-25 

1.6  REDUCED  LINEAR  COEFFICIENT  OF  THERMAL 

EXPANSION .  1-30 

TABLES 

1-1  COMPUTATION  SHEET  -  STRESS  RELAXATION 

MODULUS,  E(t) .  1-7 

1-2  COMPUTATION  SHEET  -  SHIFT  FACTOR  ax .  1-11 

1-3  COMPUTATION  SHEET  -  STRESS  RELAXATION 

MODULUS,  E(t) .  1-13 


1-1 


LOCKHKKD  PROPULSION  COMPANY 


FIGURES 


Page 

1-1  Tab  End  Tensile  Specimen,  Schematic .  1-4 

1-2  Tab  End  Tensile  Specimen  Under  Test .  1-5 

1-3  Force  and  Displacement  versus  Time,  Ramp  Stress 

Relaxation  Test .  1-6 

1-4  Stress  Relaxation  Data  for  Temperatures  Between  -110°F 

and  +174°F,  Polycarbutene-R  Propellant  Batch  7958  ...  1-8 

1-5  Temperature -Time  Shift  Relationships,  T«rT8 .  1-9 

1-6  Temperature -Time  Shift  Relationships,  T>TS .  1-10 

1-7  Master  Stress  Relaxation  Modulus  at  T  =  70°,  150°,  and 
-68®F  and  Shift  Factor,  Polycarbutene-R  Propellant 
Batch  7958  .  1-14 

1-8  JANAF  Tensile, Specimen  with  Plastic  Strain  Gage .  1-16 

1-9  Typical  Tensile  Test  Force-Time  Curves .  1-17 

l- 10  Constant  Crosshead  Rate  Properties,  Polycarbutene-R 
Propellant  Batch  7958,  JANAF  Measured  Strain 
Technique .  1-19 

1-11  Time- Temperature  Data  Conversion  Procedure .  1-20 

1-12  Constant  Load  Data .  1-21 

1-13  Master  Curve  for  Constant  Load  Strength  versus  Time  to 

Rupture .  1-23 

1-14  Bond  in  Tension  Specimen .  1-26 

1-15  Bond  Test  Specimen  Showing  Cohesive  Propellant 

Failure .  1-27 

1-16  Constant  Rate  Liner/Propellant/Substrate  Bond  Data, 

Polycarbutene-R  Propellant .  1-28 

1-17  Quartz  Tube  Dilatometer .  1-31 

1-18  Reduced  Propellant  Linear  Coefficient  of  Thermal 

Expansion  versus  Temperature .  1-32 

1-2 


LOCKHEED  PROPULSION  COMPANY 


Section  1 


ENGINEERING  PROCEDURES  FOR 
PROPELLANT  PHYSICAL  CHARACTERIZATION 


The  following  subsections  describe  methods  for  determining  those  specific 
propellant  physical  properties  used  in  performing  the  grain  structural  analy¬ 
ses  given  in  Sections  2  and  4.  Engineering  techniques  for  viscoelastic 
property,  failure  criteria  and  propellant  thermal  expansion  coefficient  deter¬ 
mination  are  described  and  examples  presented.  The  methods  of  testing  and 
data  handling  described  have  been  selected  such  that  only  conventional  labora¬ 
tory  apparatus  is  required  for  execution  of  the  tests  and  only  slide  rule  or 
desk  calculation  methods  need  be  used  in  data  reduction. 


Eachof  the  mechanical  tests  described  can  be  performed  using  a  conventional 
tensile  testing  apparatus,  such  as  the  Instron  Universal  Tester.  The  deter¬ 
mination  of  the  coefficient  of  thermal  expansion  of  propellant  requires  a  linear 
dilatometer.  Alternate  and  in  some  instances  more  economical  methods  exist 
for  conducting  the  various  tests.  The  test  procedures  described  are,  there¬ 
fore,  representative  of  the  general  method,  and  illustrative  of  their  basic 
simplicity.  Description  of  preferred  methods  for  propellant  test  specimen 
preparation  and  handling  is  included  in  Appendix  A  of  the  handbook. 


General  background  information  pertinent  to  the  tests  is  discussed  in  Section 
3  and  additional,  more  specialized  methods  of  test  and  data  handling  are  de¬ 
scribed.  The  subsection  organization  of  Sections  1  and  3  correspond  numeri¬ 
cally  and  topically  for  convenience  in  cross  referencing. 
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1.1  MASTER  STRESS  RELAXATION  MODULUS  (Reference  Section  3.1) 


Purpoee:  Determination  of  the  stress  relaxation  modulus  by  ramp  strain, 
stress  relaxation  techniques.  Isothermal  and  reduced  variable  data  reduction 
methods  are  described  and  illustrated.  The  time -temperature  correspondence 
concept,  used  in  reduced  variable  treatments,  is  emphasized  as  an  engineer¬ 
ing  method  for  obtaining  the  broad  spectrum  master  relaxation  moduli  needed 
in  grain  analysis  (Section  1.1.1). 

Test  Equipment:  Instron  Universal  Tester  with  temperature  conditioning 
chamber,  or  equivalent  apparatus. 

Test  Specimen:  Machined  tab-end  tensile  specimen  (Appendix  A)  as  shown 
schematically  in  Figure  1-1: 


Figure  1-1  Tab  End  Tensile  Specimen,  Schematic 


with  dimensions 

Propellant  length  =  lQ  (in. ) 

Cross  section  area  =  h*  (in.1)  =  Aq 
in  the  unstrained  state. 

Test  Method:  The  tab-end  specimen,  conditioned  at  a  constant  temperature 
is  clamped  in  the  jaws  of  the  tester  (Figure  1-2)  and  extended  at  a  constant 
crosshead  speed  to  a  predetermined  length  where  the  crosshead  is  stopped. 
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Figure  1-2  Tab  End  Tensile  Specimen  Under  Test 


tost 
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The  decay  of  force  at  the  extended  length  le  measured  at  time  progresses. 
With  reference  to  Figure  1-1  and  Figure  1-3,  the  following  constants  for  the 
teat  analysis  are  defined: 


TIME,  t,(SEC) 


Figure  1-3  Force  and  Displacement  versus  Time, 
Ramp  Stress  Relaxation  Test 


LI 

Strain  at  time  1 1  =  ti  *  (in. /in.) 


F(1  +  «,) 

True  stress  (for  t  >  ti)  =  . . .  (psi) 


F  (l  +«,) 

Relaxation  Modulus  =  E(t)  =  . .  . . 

A  €  ■ 

O  1 

Strain  levels  «j  between  0.01  and  0.10  are  commonly  employed  in 
this  test. 

Reduction  of  Isothermal  Data:  The  relaxation  modulus  at  a  given  temperature 
for  times  greater  than  10  ti  only  is  calculated  and  plotted  on  log-log  graph 
paper.  Table  1-1  illustrates  calculation  of  the  relaxation  modulus  at  T*70*F. 
Figure  1-4  depicts  die  data  of  Table  1-1  in  addition  to  data  for  several  other 
temperatures. 
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LOG  STRESS  RELAXATION  MODULUS  E  (t)  <*•*■> 


LOG  TIME  (SEC) 


Stress  Relaxation  Data  for  Temperature 
Polycarbutene -R  Propellant  Batch  7958 


Data  for  Temperatures  Between  -110°F  and  + 
ropellant  Batch  7958 


1.1.1  Reduced  Variable  Data  Reduction 


Determination  of  the  Time -Temperature  Shift  Factor:  An  arbitrary  reference 
temperature  TB  is  selected.  Conveniently,  the  teat  temperature  nearest  the 
ambient  laboratory  room  temperature  is  chosen. 


For  data  at  test  temperatures  T  less  than  T a,  data  such  as  those  of 
Figure  1-4,  are  shitted  to  the  left.  The  relative  shift  a^p  necessary 
to  Dring  adjacent  curves  into  superposition  is  determined. 


Figure  1-5  Temperature -Time  Shift  Relationships,  T<  Ts 


For  example,  from  Figure  1-5, 

aTi  =  |^  or  log  aTi  =  log  t2  -  log  ti 
0 

aTj  =  or  log  aTj  =  log  t2  '  -  log  ti 

The  shift  factor  aT  for  temperature  Ti  is  a,^  as  given.  The  shift  factor  for 

temperature  T2  is  a,^  x  a^  (or  log  a,^  +  log  a^).  The  process  is  repeated 

for  each  temperature;  the  shift  factor  for  data  at  the  lowest  temperature  is 
given  by  the  product  of  the  a^,n  (or  the  sum  of  the  log  a,pn). 
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For  data  at  a  test-temperature  T  greater  than  T«,  data  are  shifted  to 
the  righE  The  relative  shift  aj^necesaary  to  bring  adjacent  curves 
into  superposition  is  determined. 


Figure  1-6  Temperature-Time  Shift  Relationships,  T>  Ta 

For  example,  with  reference  to  Figure  1-6, 
aTj  =  |J  or  log  aTj  =  log  t,  -  log  t4 

aT4  =  or  lo8  aT<  =  log  ts'  -  log  t4' 

The  shift  factor  aT  for  temperature  Tj  is  as  given.  The  shift  factor  for 
temperature  T4is  aTj  x  a,^  (or  log  aTj  +  log  aT<  =  log  aT). 

The  shifts  ajn  are  shown  on  Figure  1-4  and  the  associated  calculations  for 
determination  of  log  aT  are  shown  in  Table  1-2.  The  circled  notations  in 
the  "log  aTn"  column  of  Table  1-2  identify  the  shifts  shown  on  Figure  1-4. 

Conversion  of  the  data  of  Figure  1-4  to  master  curve  form  is  accomplished 
by  performing  the  shift  aT  for  the  data  at  each  temperature.  This  corres¬ 
ponds  to  division  of  the  time  scale  for  the  data  at  each  temperature  by  the 
value  of  aj  for  that  temperature.  More  conveniently,  the  operation  is  con¬ 
ducted  logarithmically  (log  t/a^,  =  log  t  -  log  a^). 
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The  calculation  of  t/aj  is  illustrated  in  Table  1-3,  Column  3  and  the  re¬ 
duced  master  relaxation  modulus  curve  is  shown  in  Figure  1-7.  The  master 
relaxation  modulus  curve  displays  the  modulus  versus  time  behavior  of  the 
propellant  at  the  temperature  Tg  =  70°.  The  relaxation  curve  for  another 
temperature  is  obtained  by  shifting  the  master  time  scale.  The  shift  for  the 
temperature  of  interest  is  achieved  by  multiplication  along  the  time  scale  by 
the  value  of  a>p  for  the  given  temperature,  or  logarithmically,  by  adding  log 
ax  to  the  master  time  scale.  Shifted  time  scales  for  temperatures  of  150°F 
and  -68#F  are  shown  in  Figure  1-7. 

By  convention,  master  data  obtained  by  reduced  variable  techniques  is  indi¬ 
cated  by  the  notation,  for  example,  E(t)  where  t  is  equivalent  to  t/a x- 
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1.2  CONSTANT  CROSSHEAD  RATE  ULTIMATE  PROPERTIES 
(Reference  Section  3.2)1 


Purpose :  Determination  of  the  strain  at  nominal  maximum  load  and  the  true 
stress  at  nominal  maximum  load  as  measured  under  constant  crosshead  rate 
test  conditions. 

Test  Equipment:  lnstron  Universal  Tester  with  temperature  conditioning 
chamber  or  equivalent  apparatus. 

Test  Specimen:  Machined  tab-end  specimen  (Figure  1-1,  Section  1.1)  or 
machined  JANAF  "dumbbell"  with  provision  for  strain  measurement  in  the 
gage  section  (Figure  1-8  and  Appendix  A).  The  tab-end  specimen  is  preferred. 
The  common  availability  of  JANAF  specimens,  however,  requires  their  pre¬ 
sent  use.  Uniaxial  failure  data  in  this  handbook  were  obtained  in  tests  of 
JANAF  samples. 

Test  Method:  The  specimen  is  fixed  in  the  grips  of  the  tester  and  extended 
at  a  constant  tester  crosshead  speed  until  rupture  occurs.  Samples  are  pre¬ 
conditioned  and  tested  at  various  temperatures.  Although  any  crosshead 
speed  may  be  selected,  two  inches  per  minute  is  a  useful  standard  laboratory 
rate . 

Data  Reduction:  The  force  versus  time  record  of  the  lnstron  test  is  inter- 
preted  for  the  nominal  maximum  load,  Fm,  point.  Figure  1-9  shows  typical 
forms  of  the  lnstron  record  with  the  interpretation  of  Fm  indicated.  The 
strain  at  Fm  is  the  strain  at  nominal  maximum  stress  <Jm.  The  strain  t  m 
is  either  read  directly  from  the  data  of  the  strain  measuring  device  (Appendix 
A)  or  for  the  tab-end  specimen  is  given  by 

«  =  f  t 

m  t  m 
o 

where  R  is  the  lnstron  crosshead  speed  (in. /sec)  and  tm  (sec)  is  the  time  at 
which  Fm  occurs.  l0  is  the  unstrained  length  of  the  propellant  bar  (Figure  1-1, 
Section  1.1).  The  nominal  maximum  true  stress  cm  for  either  specimen  is 
given  by 


m 


F 

m 

i r 

o 


(1 


+  <  ) 
m' 


The  test  procedure  conforms  to  the  Recommended  Uniaxial  Tensile  Test 
Procedure  of  Section  4.3.5  of  the  "ICRPG  Solid  Propellant  Mechanical  Be¬ 
havior  Manual",  July,  1963,  for  JANAF  Type  Specimens . 
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Figure  1-8  JANAF  Tensile  Specimen  with  Plastic  Strain  Gage 
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Figure  1-9  Typical  Tensile  Test  Force-Time  Curves 
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where  A0  is  the  unstrained  cross  section  area  of  the  specimen.  The  associ¬ 
ated  strain  rate  <  is 


<  /t 
m  m 


where  tm  is  the  time  at  which  the  nominal  maximum  force  Fm  is  observed. 
The  nominal  value  of  «'  for  a  tester  crosshead  speed  of  two  inches  per  minute 
for  a  JfARAlF  specimen  is  0.7  in. /in. /min.  Tabular  data  for  the  tests  are 
more  conveniently  represented  graphically  as  in  Figure  1-10. 

Figure  1-10  also  shows  data  for  other  strain  rates  obtained  through  use  of  the 
temperature -time  shift  relationship  for  the  propellant.  This  data  conversion 
is  necessary  in  grain  failure  analysis,  as  discussed  in  Section  2.2.  The  pro¬ 
cedure  is  as  follows  (Reference,  Figure  1-11): 


(a)  With  curves  for  the  log  a-p  shift  factor  versus  temperature  and 
ultimate  properties  versus  temperature  for  the  propellant 
(Figures  1-7  and  1-10);  take  the  ratio  of  the  reference  (experi¬ 
mental)  strain  rate  to  the  strain  rate  of  interest,  «  ,  defining 
‘og  »Tt: 

Log,  log  'experimental 

( 

r 

(b)  At  a  given  temperature  Ti,  add  log  a  to  the  log  aT  curve  at 

T»;  lr 

(c)  Proceed  horizontally  from  the  point  aq-  +  aq-r  (point  1  in  Figure 
1-11)  left  to  the  propellant  log  a>p  curve  intercepting  it  to  define 
a  temperature  T2; 

(d)  The  ultimate  properties  at  the  new  strain  rate  are  defined  as 
equal  to  those  at  temperature  Tj  as  shown; 

(e)  The  shifted  curves  are  then  constructed  point  by  point  following 
steps  (a)  through  (d)  above. 

The  steps  given  above  are  applied  similarly  in  shifting  the  data  to  obtain 
the  curves  as  shown  in  Figure  1-10. 
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Figure  1-10  Constant  Crosshead  Rate  Properties  ,  Polycarbuter.e-R 

Propellant  Batch  7958,  JANAF  Measured  Strain  Technique 
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1.3  CONSTANT  LOAD  ULTIMATE  PROPERTIES  (Reference  Section  3.3) 


Purpose:  Determination  of  the  relaxation  of  load  sustaining  capability  of  a 
propellant  with  time  and  temperature. 

Test  Equipment:  Instron  Universal  Tester  with  load  cycling  controls  and 
temperature  conditioning  chamber,  or  equivalent  apparatus. 

Test  Specimen:  Machined  JANAF"dumbbell"or  machined  tab-end  specimen. 
Either  specimen  is  acceptable  for  this  test. 

Test  Method:  A  specimen  is  fixed  in  the  grips  of  the  tester  and  elongated  at 
a  constant  crosshead  displacement  speed  to  a  preselected  load  level,  where 
the  load  level  control  of  the  tester  has  been  set.  The  load  level  is  maintained 
until  the  specimen  either  fails  or  the  test  is  terminated.  While  any  initial 
crosshead  speed  may  be  used  in  this  test,  a  two-inch  per  minute  speed  is 
standard.  Defining  t^  as  the  time  at  which  the  preselected  load  level  is  first 
attained,  tests  are  made  at  successively  lower  load  levels  such  that  failures 
are  observed  at  times  in  excess  of  10  t|,. 

Replicate  tests  in  which  the  time  to  failure  varies  between  10  tr,  and  about 
1000  Tl  are  made,  defining  a  curve  as  in  Figure  1-T2: 


LOG  LOAD 


Figure  1-12  Constant  Load  Data 


The  procedure  is  repeated  at  various  temperatures. 
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Data  Reduction:  As  indicated  above,  the  product  data  of  the  test  are  points 
defining  a  curve  of  force  versus  time  to  failure.  The  force  readings  are  re¬ 
duced  to  obtain  stress  <rr  by 

F 

r 

%  =  7T 

o 

where  Fr  is  the  force  at  failure  (rupture)  and  Ac  is  the  unstrained  cross- 
sectional  area  of  the  specimen.  Data  points  thus  reduced  are  plotted  on 
logarithmic  coordinates  and  shifted  in  accord  with  the  shift  factor  for  the 
propellant  (Figure  1-7)  to  obtain  a  master  curve  representation.  The  shifting 
procedure  is  that  used  for  reduction  of  stress  relaxation  data  (Section  1.1). 
Typical  data  are  shown  in  the  master  curve  representation  (Figure  1-13). 


Shifted  (and  unshifted)  data  from  this  test  usually  show  relatively  large  scatter. 
Since  the  data  are  of  use  as  a  failure  criterion,  the  lower  scatter  limit  (stre  ss- 
wise)  is  identified  as  the  curve  of  interest. 
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LOG  STRESS  AT  RUPTURE  (PSD 


Master  Curve  for  Constant  Load  Strength  versus 


1.4  CONSTANT  STRAIN  ULTIMATE  PROPERTIES  (Reference  Section  3.4) 


Purpose:  Determination  of  the  constant  strain  level  below  which  the  time  to 
failure  in  uniaxial  tension  approaches  infinity. 

Test  Equipment:  Instron  Universal  Tester  with  temperature  conditioning 
chamber  or  equi vale nt . 

Test  Specimen:  Machined  tab-end  specimen  (Figure  1-1,  Section  1.1)  or 
machined  JANAF  "dumbbell"  with  provision  for  strain  measurement  in  the 
gage  section  (Reference:  Appendix  A).  The  tab-end  specimen  is  preferred. 
The  common  availability  of  JANAF  specimens,  however,  dictates  their  con¬ 
tinued  use  for  this  test. 

Test  Method:  The  specimen  is  fixed  in  the  grips  of  the  tester  and  extended  to 
a  preselected  strain  level  at  a  constant  crosshead  speed.  While  any  X-head 
speed  may  be  selected,  two  inches  per  minute  is  a  useful  laboratory  standard. 
The  sample  is  held  at  the  preselected  strain  level  and  the  time  to  failure 
observed.  The  test  is  repeated  for  various  strain  levels  to  define  a  strain 
versus  time-to-failure  (or  visible  specimen  cracking)  relationship. 

Tests  as  above  are  carried  out  at  various  temperatures. 

Data  Treatment:  Data  for  this  test  method  are  frequently  difficult  to  treat. 
It  has  been  observed  that  with  quality  PBAA,  Polysulfide,  and  Nitroplastisol 
propellants  the  strain  that  can  be  supported  almost  indefinitely  without  failure 
is  within  one  or  two  percent  (of  strain)  of  the  strain  at  nominal  maximum 
stress,  «m,  as  measured  in  the  constant  strain  rate  test  when  the  same 
straining  rate  is  used  in  both  tests.  (Reference:  Section  1.2  and  3.4.)  In 
such  instances,  the  constant  strain  test  becomes  correlative  to  the  constant 
strain  rate  test  and  extensive  definition  of  the  time-to-failure  relationships 
is  not  required.  As  a  practical  matter,  the  usual  test  observation  is  that  at 
or  above  the  strain  level  «m,  specimens  fail  within  a  few  seconds;  at  one  or 
two  percent  below,  failure  is  not  observed,  in  many  cases  after  two  or  three 
months.  In  evaluation  of  a  family  of  propellants  known  to  exhibit  this  behavior 
(on  the  basis  of  extensive  investigation),  thirty-minute  duration  constant  strain 
screening  tests  are  sufficient  to  single  out  important  propellant  deviations. 
When  such  deviations  are  observed,  (i.e.,  rapidly  decaying  constant  strain 
resistance  to  failure)  the  propellant  is  usually  poor  and  reference  to  more 
detailed  consideration  is  required  (Section  3.4). 
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1.5  PROPELLANT  -  LINER  -  SUBSTRATE  BONDING  ULTIMATE 
PROPERTIES  (Reference  Section  3.5) 


Purpose:  To  determine  the  ultimate  strength  of  the  propellant-to-liner-to- 
substrate  (insulation  or  case)  bond. 

Test  Equipment:  Instron  Universal  Tester  with  temperature  conditioning 
chamber  and  load  cycling  attachment,  or  equivalent. 

Test  Specimen:  The  bond -in-tension  specimen  is  used  (Figure  1-14).  Pro¬ 
pellant is  cast  through  the  pipe  onto  the  lined  (or  otherwise  prepared)  anvil 
surface  and  cured  as  in  the  motor  application.  Materials  used  in  preparation 
of  the  anvil  surface  are  those  used  in  the  motor  application  of  interest  and 
are  applied  in  the  same  manner  as  in  the  motor.  Conveniently,  the  Teflon 
spacer  has  an  inside  diameter  of  1.13  inches  (to  give  a  bonded  cross-sectional 
area  AQ  of  1  in.2)  and  a  height  of  3/ 16  inch. 

Test  Method  1:  The  bond  jig  assembly  is  fixed  in  the  jaws  of  the  tester  and 
the  anvil  and  pipe  sections  are  pulled  apart  at  a  crosshead  speed  one-tenth  of 
that  used  in  uniaxial  tests  of  the  propellant.  The  recommended  rates  are  0.2 
and  2  inches  per  minute.  The  bond  specimen  is  pulled  to  failure  and  the  force 
required  for  failure  observed.  The  test  is  repeated  at  various  temperatures. 

Data:  The  force  at  failure  of  the  specimen  is  recorded;  for  the 
recommended  spacer  dimension,  the  force  at  failure  is  by  definition 
the  unit  failure  stress. 

The  type,  as  to  location,  of  failure  of  the  specimen  is  recorded: 

(1)  Cohesive  in  propellant 

(2)  Cohesive  in  liner 

(3)  Cohesive  in  substrate 

(4)  Interfacial 

(4a)  Liner-propellant 
(4b)  Liner -substrate 
(4c)  Substrate -case 

When  the  failure  character  is  mixed,  the  estimated  percentages  are 
recorded  [e.g.,  60/40  (l)/(4a)]  . 

In  reduction  of  this  type  of  data  for  use  in  analysis,  a  limit  condition 
iB  used  for  routine  purposes.  Tests  in  which  over  90  percent  of  the 
failures  are  found  to  be  clearly  cohesive  in  the  propellant  (Figure 
1-15)  show  the  propellant  to  be  the  strength -limiting  component,  and 
propellant  data  as  obtained  in  the  more  precise  tests  (Sections  1.2, 

1.3)  are  used  as  the  governing  failure  criteria.  In  this  case,  the 
strength  as  measured  in  the  bond  test  agrees  reasonably  well  with 
the  uniaxial  constant  strain  rate  data  except  at  low  temperatures. 
Typical  bond  test  data  are  shown  in  Figure  1-16. 
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Figure  1-14  Bond  in  Tension  Specimen 
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Figure  1-15  Bond  Test  Specimen  Showing  Cohesive  Propellant  Failure 
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When  failure  is  dominantly  interfacial  in  the  liner,  or  in  a  substrate, 
a  materialor  process  development  problem  is  indicated.  In  general, 
suitable  liners  and  lining  techniques  yield  case  bonds  superior  in 
strength  to  propellants.  When  this  is  not  the  case,  grain  analysis 
and  resolution  of  the  associated  problems  is  decidedly  not  a  routine 
engineering  job. 

Test  Method  2:  The  bond  jig  is  fixed  in  the  jaws  of  the  tester  and  pulled  to  a 
preselected  load  level  at  a  crosshead  speed  one-tenth  of  that  used  in  uniaxial 
tensile  tests  of  the  propellant  (Section  1.2,  1.3).  The  preselected  load  level 
is  maintained  until  failure  is  observed.  The  tests  are  repeated  for  various 
temperatures  and  load  levels. 

Data:  The  treatment  of  data  for  this  test  method  are  as  above  for 
load  specification  and  failure  type  description.  In  general,  the  de¬ 
cay  of  load  sustaining  capability  for  propellant  cohesive  failure  with 
time  correlates  to  the  propellant  strength  as  measured  in  uniaxial 
tension  (Section  1.3).  It  may  be  observed  that  failure  localization 
at  an  interface  or  in  a  liner  will  occur  under  constant  load  conditions 
(Method  2)  and  not  under  constant  rate  conditions  (Method  1).  The 
analysis  of  the  data  in  any  case  proceeds  on  the  basis  of  the  limit 
considerations  delineated  in  discussion  of  Method  1. 
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1.6  REDUCED  LINEAR  COEFFICIENT  OF  THERMAL  EXPANSION 
(Reference  Section  3.  6) 


Purpose:  Determination  of  the  reduced  propellant  linear  coefficient  of  ther- 
mal  expansion. 

Test  Equipment:  Quartz  tube  dilatometer,  associated  temperature  condition- 
ing  and  measurement  apparatus,  and  electronic  recorders  (Figure  1-17),  or 
equivalent  apparatus.  Reference  reports  LPC  556-F-l,  LPC  578-F-l. 

Test  Specimen:  Cast  or  machined  propellant  cylinder  of  3/8-inch  diameter 
by  1  to  2  inche s  in  length  (/i)  or  rectangular  equivalent. 

Test  Procedure;  The  change  in  length  is  measured  with  temperature  of  the 
specimen  under  cooling  (heating)  rates  of  60°F/hour  or  less.  For  a  given 
sample  of  propellant,  cooling  and  heating  cycles  are  repeated  until  the  identi¬ 
cal  specimen  length  versus  temperature  variation  is  observed  on  successive 
heating  and  cooling  cycles.  Details  of  the  test  procedures  may  be  found  in  the 
previously  cited  references.  Samples  are  desiccated  and  maintained  in  a 
desiccated  environment  during  the  test. 

Data:  The  change  in  length  per  unit  length,  =  H/t],  is  plotted  versus 
temperature,  taking  A  L  =  0  and  -ti  at  a  reference  temperature  Ti.  Ti  is 
selected  to  be  the  zero  strain  temperature  for  a  motor  cast  with  the  subject 
propellant.  Typical  values  for  Ti  are  cure  temperature  plus  15°F  for  PBAA 
propellants  and  cure  temperature  plus  22°F  for  Nitroplastisol  propellants. 
Data  for  Polycarbutene  propellant  are  shown  in  Figure  1-18  referenced  to 
Tj  =  1 55°F. 
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Section  2 


ENGINEERING  METHODS  FOR  GRAIN  STRUCTURAL  ANALYSIS 


In  preparing  to  analyze  a  solid  propellant  rocket  grain,  it  is  necessary  as 
well  as  expedient  to  establish  the  assumptions  under  which  the  analysis  will 
be  conducted.  Notwithstanding  certain  obvious  deficiencies,  particularly  in 
connection  with  failure  theory,  the  medium  is  considered  as  isotropic,  homo¬ 
geneous,  and  continuous.  The  practical  objections  to  these  assumptions  are 
based  upon  the  fact  that  the  viscoelastic  elastomer  (consisting  of  the  order  of 
30percent  by  volume)  is  merelya  binder  in  which  various  considerably  harder 
particles  are  imbedded.  Thus,  a  homogeneous  medium  does  not  exist- -nor 
does  an  isotropic  medium.  Nevertheless,  it  is  assumed  that  on  the  average 
macroscale  there  does  exist  an  equivalent  medium  of  this  type.  For  most 
analyses  this  approximation  is  satisfactory,  although  the  assumption  can  be 
seriously  in  error  in  fracture  or  tearing,  where  the  origin  of  failure  begins 
on  the  microscale.  Next,  the  assumption  of  continuity  is  not  always  fulfilled 
because  it  implies  that  there  is  always  a  bond  between  the  various  solid  filler 
elements  and  the  elastomeric  binder.  The  pullaway  effect  is  well  established, 
wherein  excessive  tensile  stress  will  cause  the  filler-binder  adhesion  to  part. 
Nevertheless,  to  conduct  present  analyses  it  is  customary  and,  at  least  tem¬ 
porarily,  necessary  to  assume  an  isotropic,  homogeneous  continuum. 

The  second  assumption  is  that  the  strains  will  be  sufficiently  small  to  allow 
use  of  infinitesimal  deformation  theory.  For  the  loads  and  geometries  used 
in  present  motors,  strains  of  30  percent  may  be  computed  from  infinitesimal 
theory,  which  certainly  pushe s  the  limit  of  validity  for  this  assumption.  Alter¬ 
nately,  even  without  viscoelastic  effects,  finite  strain  analysis  is  far  from 
simple.  Considering  the  widespread  knowledge  of  infinitesimal  deformation 
theory  and  the  relative  ease  of  application,  it  is  considered  appropriate  to 
begin  at  this  point. 

As  a  third  assumption,  it  is  appropriate  to  neglect  inertia  forces  due  to 
straining,  during  ordinary  viscoelastic  deformations.  The  inertial  forces 
are  usually  highly  damped  and  exponentially  decaying,  although  in  certain 
cases--such  as  possible  stress  wave  propagation  through  the  grain  due  to 
burning  phenomena-- this  particular  assumption  ultimately  may  have  to  be 
re-evaluated. 

The  following  subsections  present  engineering  methods  for  grain  structural 
analysis  for  the  situations  of  primary  interest  in  preliminary  grain  design 
evaluation.  As  applicable,  the  propellant  physical  data  of  Section  1  are  used 
in  determination  of  grain  strains  or  deflections  and  grain  failure  conditions. 
The  problem  format  has  been  selected  to  simplify  rapid  desk  calculator  or 
slide -rule  calculations. 

The  analyses  presented  in  this  section  are  of  basic  importance  both  in  pre¬ 
liminary  and  final  design  evaluations.  Background  discussion  and  additional 
grain  structural  design  analysis  information  are  presented  in  Section  4. 
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2.1  THERMAL  STRAIN  ANALYSIS  -  PROBLEM  1  (Reference  Section  4.1) 


Purpose:  To  calculate  the  maximum  equilibrium  thermal  hoop  strain  «q'  8(a) 
at  the  port  in  finite  length  case -bonded  circular  and  star  port  rocket  grainB. 


Governing  Equation: 


x* =  iZ  Li 

Notes:  2.12  2.13 


1  +  v 

c 


AT) 


X2 


2.14  2.15 


2.16 


c  ,  \  Xj 

‘e(a)  =  1+  In  (T'+ Xj") 
2.11 

‘e(a)  =  «g  (a)  K. 

2.17 


‘q  (a)  =  Maximum  hoop  strain  at  the  midpoint  of  a  circular  port, 
finite  length  grain 

«g  (a)  =  Maximum  hoop  strain  at  the  midpoint  of  a  star  perforated, 
finite  length  grain 

Xj  =  Maximum  hoop  strain  at  the  midpoint  of  a  circular  port 
finite  length  grain,  infinitesimal  deformations 

X  =  b/a  =  grain  0.  D.  /i.  D. 

F  =  Parr  end  effect  correction  for  finite  length  grains 


=  Star  valley  stress  concentration  factor 

cfp  =  Propellant  reduced  linear  coefficient  of  thermal  expansion 

v  =  Case  Poisson's  ratio 
c 

Vp  =  Propellant  Poisson's  ratio,  taken  as  Vi 

a  Case  linear  coefficient  of  thermal  expansion 

AT  =  Temperature  increment  from  grain  zero  strain 
tempe  rature 
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Ce  (“)  Kj 


THERMAL  STRAIN  ANALYSIS  PROBLEM-I 


Kj  =  STAR  VALLEY  STRESS  CONCENTRATION  FACTOR 


C2  =  3/2  P,  P  =  PARR  CORRECTION  FOR 
FINITE  LENGTH  GRAINS 


PROPELLANT  r'olyCQ 

BATCH  _ Z9£2 _ 

MOTOR  _ XStA  ~  £■ 


'ffilvr.arbu  ienet  IP 


NOTE  2.16 


NOTE  2.12 


in  =  Naperian  logarithm 

Superscript  s  =  star  grain 
Superscript  c  =  circular  port  grain 

The  procedures  used  in  solving  this  equation  are  illustrated  on  the  calculation 
sheet,  Table  2- 1 .  as  annotated  below. 

NOTES: 

2.1 1  This  operation  corrects  for  the  finite  deflection  of  the  port  byan  approxi¬ 
mation-technique  that  yields  an  accuracy  of  approximately  1  percent  in  the 
calculated  strain  values  at  strains  up  to  25  percent. 

2.12  The  term 

!  p  *  c. 


includes  the  end  condition  correction  factor  F which  is  the  ratio  of  strain  in  a 
finite  length  grain  to  that  in  an  infinite  length  grain.  F  is  obtained  from  Refer¬ 
ence  Figures  D-l  and  D-2  (Appendix  D)  for  a  grain  length  L  and  diameter 
D=  2  band  of  0.  D.  tol.D.  ratio  ^  =  b/a.  In  a  grain  head  bonded  atone  end, 
obtain  F  equivalently  by  doubling  the  actual  grain  length. 

2,1  3  ffp,  the  reduced  thermal  coefficient  is  equal  to  the  change  in  length  per 
unit  length  in  the  propellant  (unconstrained)  for  the  temperature  increment 
AT  as  referenced  to  the  temperature  Ti,  the  zero  strain  temperature  of  the 
grain  (see  Note  2.15).  Values  for  Op  are  entered  directly  on  the  calculation 
sheet,  Table  2-1,  from  Figure  1-18,  Section  1.6  for  the  temperature  at  which 
strain  is  to  be  calculated. 


2.14  The  term 


1  +  v 

c 


TTv 

P 


or 

c 


1  +  V 

°c  *  C‘ 


compensates  for  the  effect  of  the  thermal  contraction  of  the  case.  For  a  steel 
case,  a  typical  value  of  Ci  is  5.05  x  10"*  in./in./°F. 

2.15  All  calculations  are  referenced  to  the  motor  zero  strain  temperature  Tj. 
The  zero  strain  temperature  is  approximately  the  temperature  at  which  the 
grain  inside  diameter  is  equal  to  the  mandrel  outside  diameter.  Typical 
values  for  Ti  are: 

PBAA,  Polycarbutene:  Cure  temperature  +  15°F  =  Ti 

Nitroplastisols:  Cure  temperature  +  22®F  =  Ti 

AT  is  the  grain  temperature  increment  from  Ti, 
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2.16  Inclusion  of  the  subtractive  term  c?p  in  the  solution  allows  prediction  of 
the  observed  strains  in  a  grain.  To  reference  the  grain  strain  calculation  to 
the  same  base  as  the  failure  data  (Section  1.2),  omit  the  subtraction  in  the 
computation  as  indicated  on  the  calculation  sheet  (  i.  e.,  the  observed  strains, 
actual  port  diameter  to  mandrel  diameter,  are  less  than  the  stress  producing 
strains  by  the  amount  mentioned  above  )  . 

2.17  The  star  valley  stress  concentration  factor  Kj  accounts  for  the  strain 
field  perturbation  at  the  star  valley  of  a  grain.  For  a  circular  port  grain 
Kj  =  1.  Kj  is  found  through  the  following  relationship: 

K  .  =  H  (iL-L_L) 

1  2X* 


i  ,/  O.D. 

\  =_  b/a  =  grain  ^  p  ■ 


from  Reference  Figures  D-3  through  D-8,  for  various  Btar  grains.  The  data 
are  from  Fourney  and  Parmerter,  1961,  with  the  permission  of  the  authors. 
Figure  D-3  specifies  the  geometric  nomenclature.  The  value  of  Kj  obtained 
from  these  figures  is  valid  only  for  star  valley  shapes  in  which  the  sides  of 
the  star  are  approximately  parallel  and  the  bottom  of  the  star  valley  is  an  arc 
of  a  circle  in  cross -section.  Complex  star  shapes  require  photoelastic  deter¬ 
mination  of  the  stress  concentration  factor. 
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2.2  THERMAL  STRAIN  FAILURE  ANALYSIS- PROBLEM  1A  (Reference 
Section  4.2) 


Purpoae:  To  predict  the  low  temperature  storage  limit  for  circular  or  star 
port  case -bonded  rocket  grains. 

Propellant  Data:  Strain  at  nominal  maximum  stress  «m  as  determined  in 
accordance' with  Section  1.2;  also,  reference  Section  1.4. 

Procedure:  Determine  the  approximate  thermal  equilibrium  time  teq  of  the 
subject  motor  (about  one  day  for  a  1-foot  motor,  about  three  days  for  a  3-foot 
diameter,  about  three  weeks  fora  10-foot  motor,  uninsulated).  Taking  the 
strain  (from  calculations  in  Problem  1)  in  the  region  of  the  lowest  temper¬ 
ature  of  interest,  define  the  mean  strain  rate  as 


c 


C,s 

9  (a) 

t 


R 


y 


and  determine  the  propellant  *m  versus  temperature  by  the  procedure  of 
Section  1.2  and  plot  on  Figure  2-1  for  the  strain  rate  Ry.  Superimpose  on 
Figure  2-1,  the  strain  versus  temperature  requirements  as  obtained  in  Prob¬ 
lem  1.  The  intercept  of  the  two  curves  defined  indicates  the  probable  crack¬ 
ing  temperature  of  the  grain  to  within  about  ±10®F  for  the  worst  condition  of 
monotonic  cooling  (i.e.,  the  fastest  rate  possible  by  normal  conditioning-box 
environments).  A  minimum  safety  factor  of  1.5  (strain  capability  4  strain 
requirements)  at  the  low  temperature  cycling  limit  of  the  grain  is  recommended 
for  motors  to  be  subjected  to  extended  storage,  cycling,  vibration,  thermal 
shock,  etc. 1  Sustained  thermal  cycling  of  a  motor  usually  will  cause  failure 
by  web  cracking  at  temperatures  higher  than  those  predicted  above  but  within 
the  limits  imposed  by  the  recommended  safety  factor  (Section  4.2). 


1  Superposition  of  pressurization  strains  and  thermal  strains  for  pressuri¬ 
zation  failure  analysis  is  discussed  in  Section  2.5. 
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STRAIN  (PERCENT) 


Figure  2-1  Thermal  Strain  Failure  Analysis,  Problem  1A 
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2.3  GRAIN  SLUMP  DEFLECTIONS  -  PROBLEM  2: 


Purpose:  Calculation  of  grain  slump  deflections  versus  time  from  acceler- 
ation  imposed  coaxial  with  the  propellant  grain.  Circular  port  grains  enclosed 
in  rigid  cases  with  free -free  and  free -fixed  end  conditions  are  considered 
(Details  9,  10.  11.  Tables  2-2  and  2-3). 

Governing  Equation: 


U 


r,  z 


=  (h 


P  ng  b2 


U  =  Deflection;  Sub  r:  radial.  Sub  z:  axial 


~Tb. 


Constant  for  solution  for  maximum  axial 
deflection  of  free -free  grain  (Figure  D-10) 


C,  Constant  for  solution  for  maximum  radial 
deflection  of  free -free  grain  (Figure  D-ll) 


D,  Constant  for  solution  for  maximum  axial 
deflection  of  free -fixed  grain  (Figure  D-12) 

F,  Constant  for  solution  for  radial  deflection 
at  free  endof  fixed-free  grain  (Figure  D-13) 


IG,  Constant  for  solution  for  radial  deflection 
at  fixedendof  fixed-free  grain  (Figure  D- 14) 


p  Propellant  density 


ng  Acceleration  in  gravities 


2b  Grain  outer  diameter 


E(t)  FVopellant  master  stress  relaxation  modulus 


Solution:  The  deflections  are  calculated  as  functions  of  time  using  the  stress 
relaxation  modulus  E(t)  appropriate  to  the  temperature  of  interest  from 
Figure  1-7,  Section  1.1.  The  calculations  for  the  various  deflections  are 
illustrated  by  Tables  2-2,  2-3  and  Figures  2-2,  2-3.  The  numerical  factors 
f>  =  B,  C,  D,  F,  G  are  scaled  from  Figures  D-10  through  D-14  for  the  motor 
I.D./O.D.  ratio  1/X  =  a/b  and  grain  length  to  diameter  L/D  ratios. 

For  grains  with  irregular  ends  (i.e. ,  non-planar),  the  values  of  the  Tr>  for  the 
grain  based  on  the  L/D  for  the  bonded  cylindrical  length  and  the  port  length 
are  determined  and  the  larger  of  the  two  values  for  p  used  in  calculation. 
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DIRECTION  OF  ACC 


663-004-REV- 1 


TABLE  2-2 


H8t 
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COMPUTATION  SHEET  -  GRAIN  SLUMP  DEFLECTIONS  -  FIXED-FREE  GRAIN  -  PROB  2B 


It* 


II  £ 


1/A  =  a/b 


ur3~^  ^ 


Detail  10 

_  . 

50“  — 


Uz2=  D  pngb2 
E(r) 

Ur3  =  F  pngb2 
E  (r) 


Ur2  =  G  pngb2 
E  (r) 


UZ2  L 


GRAIN  LENGTH  =  L  =  Zoo ' 
GRAIN  DIAMETER  2b  /  0O " 


Detail  11 


FROM  REF  FIG BLzJTL,  D  =  5.  5 p  =  PROP  DENSITY  = 

FROM  REF  FIG £-/•?.  F  =  0.109JZ.  b2  = 

FROM  REF  F\G0-rf,  G  =  j2jl  /££. _  „„  / 


(GRAVITIES) 


(LB/IN.3) 


X3  =  D  pngb2  =  5~7.  £T 
X4  =  F  pngb2  =  / •  ^** 
X5  =  G  pngb2  = 


LOG  X3  =  (l_  ^  ^ 

LOG  X4  =  _  LlAA 
LOG  Xs  =  . 


LOG  E  (r)  FROM 


/—  7 


O  Aft  I 
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DIRECTION  OF  ACCELERATION 


TABLE  2-3(cont'd 
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The  method  of  solution  assumes  a  step  function  application  of  the  constant 
acceleration  ng  at  zero  time.  For  storage  slumps  zero  time  may  be  taken 
when  the  motor  casting  hardware  is  removed.  For  launch  acceleration,  zero 
time  is  taken  as  ignition  time  and  the  previous  storage  deflections  are  added 
to  obtain  total  deflections. 


The  method  of  solution  also  assumes  that  the  reciprocal  of  the  relaxation 
modulus  E(t)  is  equal  to  the  creep  compliance  D(t).  This  is  an  engineering 
approximation  of  acceptable  accuracy  as  discussed  in  Section  3.1. 
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2.4  GRAIN  PRESSURIZATION  -  PROBLEM  3  (Reference  Section  4.4) 

Purpose:  Calculation  of  the  limit  strains  at  the  port(s)  in  circular  and  star 
port  grains  case-bonded  in  elastic  cases.  1 

Governing  Equations: 

2.4.1  Grains  with  Free  Ends  (metal  case): 

3.4  b  \lp. 

- —  K. 

hE  (3  +  V2)  1 

c 

b/a  =  grain  O.  D. /i .  D.  ratio 
grain  O.  D. 
case  thickness  (in.  ) 
case  elastic  modulus  (psi) 

Star  valley  stress  concentration  factor  (Refer¬ 
ence  Section  2.17  and  Figures  D-3  -  D-8). 

p.  =  motor  operating  pressure 

2.4.2  Grains  with  Bonded  Ends  (metal  case): 

b  Xz  (0.91)p. 

«eP  <a>  =  TT - 1  Ki  =  °-27  <3  +  x  )  «ep  (a> 

(Terminology  as  for  2.4.1) 

The  above  calculations  yield  limits  for  the  maximum  tensile  strain  at  the  port 
in  a  star  or  circularly  perforated  grain  with  the  appropriate  end  fixations. 
For  grains  with  one  fixed  and  one  free  end,  the  strains  Cg'p(a)  and  «Qp(a)  indi¬ 
cate,  respectively,  upper  and  lower  limits  for  the  maximum  strain.  Propel¬ 
lant  compressibility  is  not  taken  into  account  in  the  above  solutions,  but  is 
ordinarily  of  negligible  importance  in  moderately  designed  grains  (<65$  web 
fraction,  operating  pressure  <1000  psi).  Methods  for  calculating  strain  com¬ 
ponents  arising  from  propellant  compressibility  are  discussed  in  Section  4.4. 


\  = 

2b  = 

h  = 

E  = 
c 

K.  = 

i 


Typical  strain  values  are  illustrated  in  Table  2-4. 


1  In  conventional  steel  case  designs  of  web  fractions  up  to  60  -  80  percent, 
pressurization  stresses  are  compressive  only.  This  subject  is  discussed 
in  Section  4.4  in  greater  detail. 


2-16 


LOCKHKKO  PROPULSION  COMPANY 


S/I  CALCULATION  SHEET 


TABLE  2-4 


■■■■■■ 


■■■■■■■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■■■■ I 


IV'JL'M'  Ci 


HHHHIIMHH  bbbbe,  jr  w.  1  r.  •want' a  bbb  ■ 


■anal 
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2.5  PRESSURIZATION  STRAIN  FAILURE  ANALYSIS  -  PROBLEM  3A 


Purpose:  To  predict  the  low  temperature  firing  limit  for  circular  or  star 
port  case  bonded  grains. 

Propellant  Data:  Strain  at  nominal  maximum  stress,  €m.  as  determined  in 
accordance  with  Section  1.2. 

Procedure:  The  limit  strain  of  pertinence  as  determined  from  Section  2.4  is 
calculated.  The  time  to  full  pressurization  of  the  motor,  tp,  is  determined 
and  the  approximate  mean  pressurization  strain  rate  Rp 


R 

P 


determined.  The  propellant  strain  at  maximum  stress  appropriate  to  the 
strain  rate  Rpis  plotted  versus  temperature  in  accordance  with  the  procedures 
for  time  temperature  shifting  (Section  1.2).  The  low  temperature  firing  limit 
is  determined  as  the  temperature  at  which  the  thermal  strain  to’ 8(a)  (Section 
2.1)  and  the  pressurization  strain  C0p(a)or  <0p(a)add  viscoelastically  to  equal 
the  available  strain  capability  of  the  propellant.  The  superscripts  c  and  s 
indicate  circular  and  star  port  grains,  as  before.  The  addition  is  performed 
as  follows: 


C  •  S  /  \ 

*0  <a) 

«  no  + 

m  y' 


*eP 


i 


where  «m(Ry)is  the  strain  capability  of  the  propellant  at  the  maximum  cooling 
rate  (Reference  Section  2.2)  and  «m(Rp)  is  the  strain  capability  of  the  propel¬ 
lant  at  the  pressurization  strain  rate,  both  versus  temperature.  Illustrative 
calculations  are  presented  in  Table  2-5  and  Figure  2-4. 

It  should  be  noted  that  the  above  definition  of  the  firing  limit  assumes  a  cumu¬ 
lative  damage  criterion.  In  practice,  a  minimum  firing  safety  factor  of  1.5 
as  determined  by  the  above  procedure  has  been  found  to  be  conservative  and 
is  recommended.  The  subject  is  further  discussed  in  Section  4.5. 
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TABLE  2-5 
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Section  3 

BACKGROUND  AND  MATHEMATICAL  PROCEDURES 


The  success  of  an  engineering  structural  analysis  results  from  application  of 
sophisticated  theoretical  concepts  and  from  use  of  accurate  physical  data  for 
the  materials  of  construction.  The  physical  data  are  of  two  essential  kinds: 
data  that  define  the  relationship  between  force  and  deformation,  and  data  that 
provide  criteria  for  prediction  of  structural  failure  in  given  force -deformation 
fields.  In  direct,  but  more  restrictive  definition:  the  material  stress-strain 
law  and  failure  criteria.  As  previously  indicated,  the  stress -strain  relation¬ 
ship  for  solid  propellants  must  intrinsically  include  time  and  temperature  as 
essential  defining  parameters.  This  requirement  is  true  also  for  definition 
of  propellant  failure  criteria. 

Section  3.1  of  the  handbook  describes  the  methods  of  linear  viscoelasticity 
and  temperature -time  reduced  variables  as  specifically  applicable  to  propel¬ 
lants.  Because  of  its  relative  novelty  in  engineering  analysis,  considerable 
background  in  this  area  is  included.  Engineering  methods  for  mathematical 
viscoelastic  properties  identification  are  described  and  numerical  examples 
presented. 

Sections  3.2  through  3.6  provide  background  information  pertinent  to  the  engi¬ 
neering  techniques  described  in  Section  1.2  through  1.6.  Section  3. 7  describes 
and  discusses  the  biaxial  diametral  test  technique  which  is  useful  both  in 
obtaining  failure  criteria  information  and  in  application  as  an  economical 
primary  test  method  for  physical  screening  of  toxic  or  rare  propellant  for¬ 
mulations.  Section  3, 8describes  the  analog  test  technique  as  applied  to  strain 
evaluation  of  propellants  in  subscale  motors. 
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3.1  PROPELLANT  VISCOELASTIC  PROPERTIES 


To  predict  the  response  of  viscoelastic  materials  to  applied  stress,  it  is 
necessary  to  know  the  elastic  and  viscous  parameters  of  the  material  as  a 
function  of  time,  rate,  and  temperature.  In  principle  these  parameters  may 
all  be  combined  into  a  generalized  stress -strain  law  such  that  the  strain 
(streBs)may  be  calculated  or  deduced  for  an  applied  stress  (strain)  as  a  func¬ 
tion  of  these  parameters.  When  it  is  justified  to  assume  this  behavior  is  of 
linear  form,  based  upon  experimental  evidence,  this  association  between 
stress  and  strain  can  be  viewed  as  a  transfer  function  having  the  property 
which,  when  multiplied  by  a  linear  functional  of  strain  (stress),  generates  the 
associated  linear  functional  of  stress  (strain).  The  more  familiar  behavior 
of  metals  is  considered  first. 


3.1.1  Operational  Concept  for  Stress-Strain  Relations 


When  a  uniaxial  steel  tensile  specimen  at  a  moderate  temperature  is  stretched, 
it  is  observed  that  the  instantaneous  strain  is  very  nearly  proportional  to  the 
applied  stress  up  to  the  yield  point;  and  that  if  the  bar  is  unloaded  the  same 
law  is  followed.  Consequently,  the  material  obeys  Hooke's  law  in  tension  up 
to  the  yield  point  and 

<r  =  Ee  (3-1) 

in  which  stress  a  ,  and  strain  e,  are  based  on  the  original  specimen  dimen¬ 
sions  as  in  the  usual  engineering  sense.  The  proportionality  constant,  E, 
which  represents  the  slope  of  the  stress-strain  curve,  is  called  Young's 
modulus  or  the  tensile  modulus .  This  law  can  be  written  also  in  the  inverse 
form 


Dcr  =  e  (3-2) 

The  proportionality  constant  D  is  then  commonly  referred  to  as  the  tensile 
compliance.  It  has  been  found  experimentally  that  the  above  law  holds  for 
many  materials,  particularly  metals,  as  long  as  the  strains  are  small. 

It  is  evident  that  if  a  material  which  obeys  Hooke's  law  is  held  at  a  constant 
strain,  the  stress  also  remains  unchanged  with  time.  However,  when  a 
viscoelastic  tensile  specimen  at  room  temperature  is  stretched  and  held  at  a 
constant  strain  eQ  (stress  relaxation  test),  the  stress  <r  (t)  necessary  to  main¬ 
tain  this  elongation  decays  with  time.  In  other  words,  the  tensile  relaxation 
modulus,  E(t)  =  tr(t)/«0,  decreases.  This  situation  is  illustrated  in  Figure 
3- 1(a), behavior  at  constant  stress  is  shown  in  Figure  3- 1(b). 

In  addition  to  strong  time  dependence,  the  mechanical  properties  are  greatly 
affected  by  temperature.  Below  a  temperature,  Tg.  defined  as  the  glass 
transition  temperature,  the  propellant  is  essentially  glassy  and  behaves  as  a 
brittle  body  obeying  Hooke's  law.  Above  this  temperature,  however,  the 
response  is  time  dependent  and  varies  considerably  with  temperature.  This 
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Figure  3-1  Typical  Propellent  Behavior  (a)  Relaxation  Test,  (b)  Creep  Test 


<(t) 


Figure  3-2  Creep  Compliance  for  Stress -Strain  Law:  c  (t)=  [Ev  +  tivjj- ]e  (t) 
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behavior  leads  to  formulation  of  a  general  functional  relationship  between 
tensile  stress  and  strain  which  includes  both  time  and  temperature  dependence 
for  temperatures  greater  than  T„.  However,  it  has  been  found  for  many 
polymers,  particularly  plastics  ana  rubbers,  that  these  two  variables  can  be 
considered  separately  if  the  temperature  range  is  not  too  great.  Forexample, 
if  certain  material  constants  are  known  at  one  temperature,  it  is  possible  to 
predict  behavior  at  another  temperature  by  simply  shifting  the  time  scale. 
Since  composite  propellants  are  filled  rubbers,  it  is  expected  that  the  same 
rule  should  hold.  Therefore,  on  this  basis  only  the  time -dependent  behavior 
at  a  fixed  temperature  will  be  discussed,  but  will  be  followed  in  later  sections 
by  an  explanation  of  the  method  used  to  change  to  a  different  temperature. 

Expanding  now  the  simple  stress -strain  relation  given  by  (3-1)  for  a  simple 
elastic  behavior  to  a  more  general  time  dependent  behavior  one  can  write 

O,  [  or  (t)]  =  02  [  «(t)]  (3-3) 

where  Oi  and  O2  represent  algebraic  and  differential  operations  on  <r  (t)  and 
e(t).  For  example,  when  Hooke's  law  (3-1)  applies,  the  operators  are  the 
simple  constants  Oi  =  1  and  O2  =  E.  It  is  important  to  note,  however,  that 
these  operators  are  not  always  linear.  Indeed,  for  large  deformations  of 
some  metals,  a  more  realistic  elastic  is,  for  example 

<r  =  E  (1-2  ve)  In  (1  +  e)  (3-4) 

where  v  is  Poisson's  ratio.  In  this  instance,  0*  is  a  non-linear  algebraic 
operator. 

A  simple  time -dependent  extension  of  Hooke's  law  is  to  consider  stress  pro¬ 
portional  to  both  strain  and  strain  rate.  For  this  case,  (3-3)  becomes 


"■(0  =  \  +  Ev  (3_5) 

in  which  Ev  and  Tly  are  proportionality  constants.  An  important  implied 
property  of  the  differential  operator  is  that  it  is  linear  and  therefore  obeys 
many  of  the  ordinary  rules  of  algebra,  such  as  association,  commutation, 
and  superposition.  This  allows  factoring  of  (3-5)  to  yield 

*(t)  =  [  \  ar  +  evi  «(t)  (3-6) 

and  hence  identify 

O,  =  \  £  +  Ey  (3-7) 

If  a  creep  test  is  considered  in  which  a  constant  uniaxial  stress,  <r0,  is  applied 
to  a  material  following  (3-5),  the  resulting  axial  strain  can  be  calculated 
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3.1.2  Linear  Viscoelastic  Representation 


A  special  generalization  of  (3-3)  is  now  defined  such  that  Oi  and  0*  are  taken 
to  be  linear  differential  operators.  In  the  literature  such  a  stress-strain 
law  is  called  a  linear  viscoelastic  representation*,  and  for  a  simple  tensile 
test  is  written 


[  a  1-  + 
n  btn 


a‘  W  +  aoJ  CT(t)=  Cbn  “n  + . +  bl  'St  +  bo]  *  (t)  {3* 1 2) 

ot 


or  more  compactly 

[P]  or(t)  =  [Q]  e(t)  (3-13) 

where  P  and  Q  stand  for  the  bracketed  operators  Oi  and  O2  respectively,  j- 
is  a  linear  operator  that  represents  the  i_th  derivative  with  respect  to  time, 
and  a^  and  bi  are  experimental  material  constants  which  are  obtained  by 
methods  to  be  discussed  subsequently. 

The  relation  (3-12)  has  been  verified  experimentally  for  small  strains  over  a 
wide  temperature  range  for  many  composite  propellants  at  useful  engineering 
strain  levels  (LPC-556-F-1). 

Before  discussing  various  specializations  of  (3-12),  it  should  be  remarked 
that  the  same  form  of  a  stress-strain  law  is  found  to  hold  for  hydrostatic 
pressure  and  volume  change,  and  for  shear  stress  and  shear  strain  in  a 
simple  shear  test.  Thus,  the  response  of  an  element  subjected  to  hydrostatic 
pressurization  is  represented  by 

[P'J  o At)  =  [Q'l  AV&  (3-U) 

H  o 

where  £V(t)/V0  is  the  volume  change  per  unit  undeformed  volume  due  to  the 
hydrostatic  pressure  °p(t),  Similarly,  write 

(  P  "]  T(t)  =  [  Q  "]  y  (t)  (3-15) 


*  If  a  time  dependent  stress  produces  an  associated  time  dependent  strain, 
then  if  doubling  the  magnitude  of  the  stress  holding  the  mode  shape  of  the 
time  variation  the  same  also  doubles  the  strain  magnitude  without  changing 
the  shape  of  its  time  dependence,  the  material  is  said  to  be  linearly 
viscoelastic. 
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in, which  T(t)  is  a  shear  stress  and  y(t)  the  corresponding  shear  strain.  P  , 
Q  **,  P  *  Q  are  of  the  general  form  of  P  and  Q  shown  in (3-12),  but,  of  course, 
with  different  experimental  material  constants.  Equations  (3-13),  (3-14),  and 
(3-15)  are  analogous  to  the  elastic  stress-strain  laws,  since  for  an  elastic 
body  undergoing  small  deformations 


simple  uniaxial  tension: 

o  =  Ec  ; 

(E  =  Young's  modulus) 

hydrostatic  pressurization: 

_  .  K  *v 

p  T 

;  (K  =  bulk  modulus)  (3-16) 

shear: 

T  =  p  y; 

(p  =  shear  modulus) 

Similarly  for  a  linear  viscoelastic  material 

simple  uniaxial  tension: 

<j(t)  =  g  ■ 

i  (t) 

hydrostatic  pressurization: 

a  (t)  = 

P  P  ' 

^ r* ~  (3-17) 

vo 

shear: 

T(t)  =  Ql 

p" 

Y(t) 

where  the  association  of  E  with 

Q/  P,  etc. , 

is  obvious. 

It  maybe  remarked  that  a  purely  mathematical  approach  to  linear  theory  does 
not  restrict  the  form  of  the  P,  Q  operators.  However,  Biot  (1954)  has  used 
irreversible  thermodynamics,  to  show  that  the  coefficients  are  restricted  in 
such  a  way  that  the  ratios  Q  / P  ,  Q  /P  ,  and  Q/P  must  be  identical  with  the 
transfer  functions  for  mechanical  models  consisting  of  springs  and  dashpots*  . 
Since  the  operator  ratios  in  (3-17)  are  restricted  to  be  of  a  definite  form,  it 
is  often  convenient  to  retain  this  form  rather  than  multiplying  all  terms  out 
and  separating,  P  and  Q.  Therefore,  E(t),  K(t)  and  p  (t)  will  be  written  in 
place  of  Q  */P  *  Q^P*.  Q/Pand  thus  maintain  this  analogy  between  the  visco¬ 
elastic  and  elastic  problems  in  the  notation. 

When  it  comes  to  actually  computing  the  stresses  and  displacements  in  a 
linear  viscoelastic  body,  use  will  be  made  of  this  analogy  by  working  out  an 
associated  elastic  solution  in  terms  of  the  associated  operator  form  of  the 
material  representation. 

It  can  be  shown  that  when  a  three-dimensional  elastic  body  is  isotropic  and 
homogeneous  there  can  exist  no  more  than  two  independent  elastic  constants 
(Sokolnikoff,  1956).  In  their  natural  form,  they  are  the  bulk  modulus,  K, 
governing  the  stress-strain  proportionality  for  (hydrostatic)  dilatation  acting 

*  Mathematical  analogies  for  spring  and  dashpot  treatments  are  included  in, 
for  example,  LPC  556-F-l,  LPC  578-F-l. 
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alone,  and  the  shear  modulus,  )x  ,  governing  distortion  alone.  Similarly,  as 
a  result  of  geometric  symmetry,  only  two  independent  operator  relations  can 
exist  for  an  isotropic,  homogeneous,  linear  viscoelastic  material.  Specifi¬ 
cally,  from  (3- 16)  and  (3- 17) 

P - »p(t)  (3-18) 


K - »K(t)  (3-19) 

Engineering  analysts  frequently  find  it  convenient  to  solve  elastic  problems 
in  terms  of  Young's  modulus  and  Poisson's  ratio  which  are  related  to  K  and 
p  .  and  the  viscoelastic  associations  by 


E 


v 


9  p  K  9p  (t)  K(t) 

3R  +  p - •*  ITO  +  pTtT  = 


3K-  2p 

_  3K(t)-2p(t) 

6K  +  2p 

6K(t)  +  2(ji  (t) 

E(t) 

v(t) 


(3-20) 

(3-21) 


In  rubber  elasticity  it  has  been  conventional  to  note  that  the  ratio  of  the  bulk 
and  shear  moduli  is  very  large,  leading  to  the  simplification 

E  =  3p,E(t)  St  3p(t)  (3-22) 

u  =  1/2,  v(t)  =  1/2  ■  (3-23) 

In  the  class  of  problems  to  be  treated  in  this  handbook,  it  has  been  found  that 
the  above  simplifications  are  adequate  for  engineering  treatment  of  well  con¬ 
solidated  solid  propellants  in  moderately  designed  solid  propellant  motors. 
In  experimental  determination  of  propellant  properties,  hence,  the  measure¬ 
ment  of  E(t)  or  p  (t)  is  the  present  problem. 


3.1.3  The  Stress  Relaxation  Modulus 


The  viscoelastic  stress  relaxation  modulus  E(t)defines  a  relationship  between 
the  time -dependent  stress  o(t)  resulting  from  an  instantaneous  straining  to  a 
constant  strain  level  «Q  and  the  viscoelastic  properties  of  the  material: 


o 


Numerically,  the  relaxation  modulus  is  obtained  in  reduction  of  data  in  which 
stress  decay  after  imposition  of  strain  on  a  sample  is  observed.  Experi¬ 
mental  techniques  for  accomplishing  the  tests  and  reduction  of  the  data  there¬ 
from  are  described  in  Section  1.1  and  will  not  be  given  in  detail  here.  The 
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objective  in  the  following  example!  will  be  to  illustrate  the  resolution  of  the 
numerical  data  to  mathematical  form  in  terms  of  the  propellant  stress-strain 
law.  This  resolution  of  data  to  generalized  form  is  an  essential  step  in  pro¬ 
viding  the  mathematical  specification  of  material  properties  required  in 
solution  of  many  stress  analysis  problems.  Two  forms  of  the  stress-strain 
law  will  be  treated:  the  linear  differential  operator  form 


P  <T(t) 
P 


Q«(t) 


n 

Q  =  E 
i=o 


(3-24) 


where  « (t)  is  time -dependent  strain,  and  the  coefficients  a.,  b.  are  the  gener¬ 
alized  material  properties;  and  the  integral  form 


<r(t)  =  I  E  +  J*  •  “7-  •  -J^-l  «  (t)  (3-25) 

L  0  -5T  +  T  J 

where  the  function  H  (t)  and  Ee  are  the  generalized  material  properties. 
Without  treating  the  methods  of  solution  of  these  equations  in  detail,  noting 
that  this  information  is  given  in  LPC  556-F-l,  LPC  578-F-l,  and  GALCIT 
SM-61-5,  the  solutions  for  the  relaxation  modulus  are,  for  the  operational 
specification, 

E(t)  =  =  A  +  A,  e’®11  +  A,e"®zt  +>--+A  e'V  (3-26) 

t  o  n 

o 

and,  in  the  case  of  the  modified  power  law  approximation,  for  the  integral 
specification 

E(t)  =  =  Ceo+C,  t‘n. 

o 

t  >  tx  (3-27) 

The  immediate  objective  is,  then,  to  obtain  numerical  values  for  the  con¬ 
stants  A^,  a. ,  or  Ci,  Ca,  and  n  from  test  data. 


(a)  The  Differential  Operator  Solution 

Assume  a  step-function  uniaxial  strain  applied  to  a  viscoelastic  bar  of  pro¬ 
pellant.  Mathematically, 
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« (t) 


t  =  o 
t  >  0 


and  as  shown  in  Figure  3-3.  For  times  greater  than  zero,  stress  decay 
shown  in  Figure  3-3  may  be  specified  by  the  following  Prony  series,  which 
is  a  solution  of  the  differential  operator  equation  (Equation  3-24) :  . 

=  A  +  A,  e'"'1  +  Aj  e'“jt  +  ...+A  .  (3-28) 

«  o  n  '  ' 

o 


Proceeding  from  experimental  relaxation  data  to  calculation  of  the  coefficients 
of  a  series  of  decaying  exponential  terms  can  be  accomplished  by  graphical 
term  stripping  techniques  (Tobolsky,  1956).  The  technique  is  based  upon 
the  observation  that  after  some  undetermined  time  all  terms  in  the  series, 
except  those  with  the  smallest  coefficients,  have  become  insignificant. 
Once  this  time  has  been  determined,  the  low-order  term  may  be  determined 
and  subtracted  from  the  data.  The  procedure  is  repeated  until  all  the  terms 
have  been  extracted. 


The  collocation  method  (Schapery,  1962;  Brisbane,  19  62)  may  be  also  uti¬ 
lized.  In  this  method  the  constants  are  assumed  and  the  corresponding 
set  of  constants  A^  are  found  from  a  set  of  simultaneous  linear  equations. 
If  the  constants  are  properly  chosen,  the  resulting  set  of  linear  equations 
can  be  solved  readily  by  desk  calculations. 


An  example  illustrates  the  collocation  method.  Data  points  taken  from  a 
typical  stress-relaxation  curve  (Figure  3-4)  at  one  decade  time  intervals  are: 


Time  (seconds) 


Modulus  (psi) 


1 

10 

100 

1000 


635 

475 

345 

265 


"n"  of  Equation  (3-28)  is  set  equal  to  four  (one  per  decade  of  time  plus  l)and 
the  value  of  AQ  is  found  to  be  200  psi.  Equation  (3-28)  now  reduces  to 


ill 


A.o 

l 


or-.  t 
i 


-  A 


=  £&  . 
€ 


200. 


(3-29) 


If  the  are  now  arbitrarily  chosen  such  that 
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3  Stress -Strain  Relationships  for  Step  Function 
Strain,  Stress  Relaxation  Test 
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STRESS  RELAXATION  MODULUS  E  <t> 


Figure  3-4  Stress  Relaxation  Data 
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where  is  the  time  at  each  data  point,  then 

O'!  =  0.5 
Oz  =  0.05 
0's  =  0.005 
o4  =  0.0005 

Using  these  values.  Equation  (3-29)  may  be  evaluated  for  the  coefficients  A^. 
The  matrix  representation  of  the  simultaneous  linear  equations  obtained  by 
writing  Equation  (3-29)  for  each  data  point  is 


Evaluating  the  exponentials. 


-0.005 

e 

py 

1.000 

-0.005 

e 

fy 

1.000 

-0.05 

e 

py 

0.951 

-0.5 

e 

py 

0.607 

-5 

e 

py 

0 

-50 

e 

py 

0 

-500  n 

e  ea  0 

the  matrix  reduced  to 
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Solving  the  resultant  set  of  simultaneous  linear  equations  gives  the  following 

values  for  A. 

1 


A,  = 

165 

A2  = 

165 

Aj  = 

71 

A*  = 

107 

and  the  Prony  series  is 

E(t)  =  200  +  165e‘0,5t  +  165e"0,05t  +  71 e'0,005t  +  107e‘°'0005t  (3-30) 


which  fits  the  experimental  data  with  an  accuracy  of  3  percent.  The  accuracy 
of  the  fit  can  be  improved  by  taking  the  arbitrary  relaxation  coefficients  at 
smaller  intervals.  In  studies  of  propellant  data,  one  exponential  per  decade 
of  time  has  been  found  to  yield  sufficiently  accurate  results  for  grain  stress 
analysis  calculations. 

Thus,  the  numerical  representation  of  the  stress  relaxation  modulus  is  ob¬ 
tained.  The  above  procedure  is  generalized  for  consideration  of  more  exten¬ 
sive  data  time  wise  by  increasing  the  number  of  terms  of  the  type 


A.  e 
x 


t 


in  the  Prony  series.  A  complete  material  specification  usually  requires 
about  15  to  20  such  terms  and  is,  on  the  practical  level,  a  one-hour  desk 
calculator  job.  Generalization  to  the  operator  equation  form  requires  only 
algebraic  manipulation  of  the  coefficients  of  Equation  (3-30);  the  relationships 
are  shown  in  Appendix  B. 


(b)  The  Integral  Modified  Power  Law  Solution 


Assume  a  step  function  uniaxial  strain  is  applied  to  a  viscoelastic  bar  of  pro¬ 
pellant.  It  is  desired  to  determine  the  coefficients  CB,  Ci  and  n  of 

E(t)  =  =  C.  +  Cx  t'n  (3-31) 

o 

the  constant  C,,  the  value  of  the  modulus  at  t  =  “is,  in  this  instance,  zero. 
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Replotting  the  data  of  Figure  3-4  on  logarithmic  coordinates  yields  a  straight 
line  function  (Figure  3-5).  The  coefficient  Ci  is  the  intercept  value  at  t  =  1 
and  the  coefficient  n  is  the  arithmetic  slope  of  the  curve,  yielding: 


E(t)  =  =  635  t‘0,170 

(3-32) 

1000  >t  >  1  • 

It  should  be  noted  that  the  modified  power  law  form  used  here  for  illustration 
is  valid  only  within  the  time  range  where  the  plot  of  log  E(t)  versus  log  t  is 
well  approximated  by  a  straight  line. 


(c)  General  Considerations 


The  preceding  examples  serve  to  illustrate  the  conversion  of  numerical 
viscoelastic  data  to  useful  forms  of  the  isothermal  stress-strain  law  for  the 
simplest  of  conceptual  test  circumstances,  i.e.,  the  step  function  strain  in¬ 
put  to  a  stress  relaxation  test.  In  experimental  practice,  the  step  function 
strain  input  is  neither  possible  nor  necessary.  Perhaps  the  present  simplest 
method  of  viscoelastic  characterization  is  the  ramp  strain,  stress  relaxation 
test,  which  can  be  easily  and  accurately  executed  using  conventional  tensile 
testers  such  as  the  Instron.  In  the  ramp  test,  strain  is  applied  to  a  tensile 
specimen  at  a  constant  rate  until  a  predetermined  strain  level  is  obtained,  at 
which  straining  is  stopped  and  stress  decay  with  time  observed.  Rigorously, 
the  appropriate  solutions  for  the  ramp  relaxation  curve  must  take  into  account 
the  stress-strain  history  during  the  ramp  portion  of  the  test  for  either  of  the 
two  forms  of  the  stress-strain  law  discussed.  On  a  practical  basis,  however, 
the  viscoelastic  memory  of  the  propellant  can  be  regarded  as  a  fading  one  in 
which  the  initial  ramp  strain  condition  is  forgotten  by  a  time  approximately 
10  times  that  of  the  ramp  duration.  Defining  the  ramp  duration  time  to  be 
ti,  the  interpretation  of  the  relaxation  curve  for  t>10tj  is,  for  engineering 
purposes,  identical  to  that  for  the  step  function  strain  input  as  discussed. 
While  there  is  loss  of  mathematical  rigor  in  this  interpretation,  numerical 
inaccuracy  is  less  than  one  percent.  The  preceding  example  may  be  inter¬ 
preted,  hence,  as  applying  either  to  a  step  function  strain  input  resolved 
from  time  zero  or  a  ramp  strain  input  resolved  for  t>10ti. 

Treatment  of  viscoelastic  data  for  wider  time  ranges,  such  as  are  obtained 
from  reduced  variable  representations,  proceeds  using  the  methods  previ¬ 
ously  illustrated.  For  the  data  shown  in  Figure  1-7.  Section  1.1,  the  Prony 
Series  coefficients  for  the  70°F  time  scale  are  shown  in  -Table  3-1.  The 
modified  power  law  representation  is 


E(t)  =  100  +  616 


(3-33) 
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Table  3-1 


PRONY  SERIES  COEFFICIENTS 
MASTER  RELAXATION  MODULUS  AT  70°F 


n 

E(t)  =  .7! 

'  '  1=0 

A,."*!* 

i 

a . 

1 

A 

i 

0 

0 

1.  00 

X 

102 

1 

5 

X 

10‘» 

1.98 

X 

101 

2 

5 

X 

10‘7 

1.  19 

X 

10* 

3 

5 

X 

10‘6 

3.95 

X 

10> 

4 

5 

X 

10‘5 

4.25 

X 

101 

5 

5 

X 

10'4 

6.50 

X 

10* 

6 

5 

X 

10'3 

9.  15 

X 

101 

7 

5 

X 

io'z 

1.69 

X 

102 

8 

5 

X 

10'* 

2.47 

X 

102 

9 

5 

X 

l 

5.71 

X 

102 

10 

5 

X 

10* 

8.  11 

X 

102 

11 

5 

X 

102 

1. 36 

X 

103 

12 

5 

X 

103 

2.56 

X 

103 

13 

5 

X 

104 

4,87 

X 

103 

14 

5 

X 

105 

8.55 

X 

103 

15 

5 

X 

10* 

1.40 

X 

104 

16 

5 

X 

107 

2.  12 

X 

104 
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3.1.4  Temperature-Time  Reduced  Variable  Treatment 


The  essence  of  the  time -temperature  reduced  variable  concept  is,  for  example, 
that  stress  relaxation  with  time  t  at  a  given  temperature  Ts  is  equivalent  to 
the  stress  relaxation  at  a  temperature  T  when  measured  on  a  reduced  time 
scale  t/a<p.  a-p  is  the  temperature -time  shift  factor.  Mathematically,  the 
relationships  in  terms  of  the  two  forms  of  the  relaxation  modulus  equation 
are 

(a)  Differential  Operator 


E(t)  =  . £  A.  e‘orit/aT 
'  1=0  1 

(b)  Modified  Power  Law 


E(t)  =  C„  +  C,  (t/aT)'n 


(3-34) 


(3-35) 


where  the  coefficients  A^,  Oj,  C®,  Ci  and  n  are  those  determined  at  the 
reference  temperature  Ts  (aj  =  1  at  Ts).  Thus  it  is  seen  that  changes  in 
temperature  correspond  physically  and  mathematically  to  changes  in  the  time 
scale  of  observation. 

The  reduced  variable  treatment  used  in  this  handbook  differs  somewhat  from 
the  WLF  method  (Williams,  Landel  and  Ferry,  1955)  in  that  an  absolute 
temperature  dependence  of  the  relaxation  modulus  is  ignored.  Discussions  of 
this  pointare  available  in  the  above  cited  reference  as  wellas  inLPC  556-F-l 
and  LPC  578-F-l.  In  the  problems  considered  in  this  handbook,  the  neglect 
of  an  absolute  temperature  dependence  of  the  relaxation  modulus  results  in 
potential  errors  in  grain  stress  analysis  that  are  trivial. 


3.1.5  Creep  Compliance 


The  stress  relaxation  modulus 


is  obtained  from  solution  of  the  operation  Equation  (3-13) 
Pa(t)  =  Q  «  (t)  or  o(t)  =  p  «  (t) 
for  the  strain  input  condition 
«  (t)  =  «o  u(t), 


(3-36) 


(3-37) 


(3-38) 
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i.e.,  a  step  strain  input  at  t  =  0.  The  analogous  solution  for  a  step  stress 
input,  solving  for  strain  is 


a(t)  =  o0u(t) 


(3-39) 


yielding  the  expression  for  the  creep  compliance  D(t) 

D(t)  =  iM  (3-40) 

which  for  the  operational  solution  is 


n 


D(t)  =  .E 
'  1=0 


B.  e'pi‘ 
1 


(3-41) 


The  Prony  Series  coefficients  Bj,  of  Equation  (3-41)  are  explicitly  related 
to  the  coefficients  Aj,  a^,  of  Equation  (3-28)  through  the  operators  P  and  Q. 
Interconversion  of  these  coefficients  maybe  done  using  the  algebraic  relation¬ 
ships  presented  in  Appendix  B;  the  subject  will  not  be  further  pursued  in  this 
section. 

It  is  important  to  notice  that  while  E(t)  is  not  exactly  the  reciprocal  of  D(t), 
to  a  satisfactory  degree  of  approximation  E(t)  can  be  regarded  as  the  re¬ 
ciprocal  of  D(t).  Hence,  the  creep  compliance  for  the  Polycarbutene-R 
propellant  used  as  an  example  in  this  handbook  may  be  numerically  estimated 
closely  enough  for  engineering  calculations  by  taking  the  reciprocal  of  the 
stress  relaxation  curve  of  Figure  1-7,  Section  1.  The  approximation  has  been 
used  in  the  solution  of  grain  slumj?  problems  in  Section  2  of  this  handbook. 

The  creep  compliance  test  is  a  useful  adjunct  in  propellant  viscoelastic  char¬ 
acterization  because  of  the  simplicity  with  which  the  test  may  be  executed 
for  extremely  long-time  test  durations.  A  test  procedure  for  creep  compli¬ 
ance  measurement  has  been  included  as  Appendix  C  to  this  handbook. 
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3.2  CONSTANT  CROSSHEAD  RATE  FAILURE  PROPERTIES 


Constant  crosshead  rate  tensile  properties  of  a  propellant  are  of  fundamental 
importance  in  both  propellant  evaluation  and  motor  structural  analysis.  The 
test  technique  is  economical  and  rapid  in  execution.  With  some  reservations, 
constant  rate  properties  provide  basic  failure  data  for  propellants  that  permit 
rapid  evaluation  of  new  formulations.  The  constant  rate  properties  are  usu¬ 
ally  the  first  available  to  the  motor  stress  analyst  when  new  formulations  are 
under  consideration  for  motor  applications.  Discussion  of  the  technical 
aspects  of  the  test  method  from  the  point  of  view  of  the  grain  structural  ana¬ 
lyst  is  therefore  important. 

The  strain  at  nominal  maximum  stress  «m  as  interpreted  from  Figure  1-9, 
Section  1.2,  provides  a  reasonably  reliable  measure  of  the  useful  strain  capa¬ 
bility  of  a  propellant  in  motor  applications  (Jones,  1963).  The  physical 
phenomenon  most  importantly  associated  with  the  strain  at  nominal  maximum 
in  composite  propellants  is  the  inception  of  extensive  oxidizer -binder  dewet¬ 
ting,  which  produces  the  change  in  properties  responsible  for  the  abrupt 
decrease  in  modulus  near 

m 

Extensive  motor  tests  have  shown  that  typical  PBAA,  PBAN,  and  Nitroplasti- 
sol  propellants  fail  when  hoop  tensile  strains  in  the  motors  exceed  the  tm  as 
measured  in  uniaxial  tension,  regardless  of  apparent  additional  elongation  to 
rupture  uniaxial  tensile  capability.  Hence,  the  implication  is  drawn  that  while 
inception  of  dewetting  may  not  cause  immediate  rupture  under  constant  rate 
uniaxial  tensile  elongation,  it  does  sounder  the  multiaxial conditions  and  slow 
rates  typical  of  thermally  cooled  motors. 

The  apparently  arbitrary  interpretation  of  the  value  of  *m  indicated  in  Figure 
1-9,  Section  1.2,  has  been  developed  by  collation  of  various  types  of  data. 
To  cite  a  few,  the  fm  as  interpreted  from  Figure  1-9  are  found  to  agree  with: 

•  The  long  term  constant  strain  capability  (LPC  556-F-l) 

•  The  strain  at  failure  under  long  duration  constant  stress  (or 
load)  conditions  (LPC  556-F-l) 

•  The  upper  strain  limit  of  useful  viscoelastic  linearity  under 
short  duration  stress  relaxation  or  creep  conditions  (LPC  578- 
F-l) 

•  The  hoop  strain  at  failure  in  scale  and  full-scale  motor  tests 
(LPC  578-F-l  ;  Hart,  1963;  Jones,  1963) 

The  latter  item  is,  of  course,  the  most  vital  of  the  types  of  correlative  infor¬ 
mation.  It  should  be  emphasized  that  these  correlations  are  valid  only  for 
comparable  conditions  of  rate  and  temperature.  The  methods  of  data  handling 
and  rate -temperature  correlation  are  illustrated  in  Sections  1.1,  1.2  and  2.1. 
Before  proceeding  to  discussion  of  the  test  method  itself,  it  should  be  noted 
that  interpretation  of  all  test  records  for  «m  is  not  possible  using  the  rules 
indicated  in  Figure  1-9.  Test  data  of  this  type  require  auxiliary  investigation, 
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moat  conveniently  by  the  constant  strain  or  constant  load  techniques  discussed 
in  Sections  3.3  and  3.4  or  in  analog  motor  tests  (Section  3.9). 

With  regard  to  the  mechanics  of  the  test  itself,  it  is  of  vital  importance  that 
a  precise  test  specimen  be  used.  The  only  generally  acceptable  method  for 
producing  sufficiently  precise  tensile  specimens  is  by  machining,  as  dis¬ 
cussed  in  Appendix  A.  Second,  an  accurate  method  of  strain  measurement 
must  be  used  in  tests  of  JANAF  tensile  specimens.  It  has  been  found  that 
the  effective  gage  length  of  a  JANAF  specimen  varies  with  propellant  type, 
temperature,  strain  rate  and  strain  level  (e.g.,  LPC-531-F,  LPC  556-F-l, 
Duerr,  1961).  Either  photographic  or  extensometer  (Appendix  A)  techniques 
can  be  used  to  measure  strain  in  JANAF  specimens  within  1  percentof  strain, 
the  extensometer  method  being  preferred  at  LPC  because  of  its  ease  and 
economy. 

In  execution  of  the  test  it  has  been  found  to  be  important  to  follow  rigorous 
specimen  storage  and  handling  precautions.  Following  completion  of  cure, 
tensile  test  specimens  are  exposed  to  undesiccated  conditions  for  only  the 
time  necessary  to  machine  them  (about  10  to  25  minutes  at  R.H.  <50  percent), 
and  during  the  one  minute  or  so  required  in  actual  execution  of  the  test  in  the 
Instron. 


3.3  CONSTANT  LOAD  PROPERTIES 


Constant  load  propellant  strength  data  provide  a  measure  of  the  long-term 
stress  capability  of  a  propellant  which  iB  useful,  for  example,  in  failure 
analysis  for  gravitational  storage  load  conditions  (Section  4-6).  This  type  of 
data  is  also  diagnostic  for  time  effects  in  failure  sometimes  not  observed  in 
the  constant  rate  tests.  For  example,  constant  load  tests  may  show  a  gradu¬ 
al  strength  decay  to  near  zero  after  an  induction  period  of  a  few  hours  in  high 
temperature  tests  that  may  not  be  evident  in  constant  rate  tests  under  the 
same  temperature  and  time  conditions.  Such  apparently  contradictory  data 
indicate  nonlinear  effects  in  a  propellant  arising,  for  example,  from  chemical 
bond  scission  in  the  binder  structure  under  sustained  load  conditions.  The 
phenomenon,  which  can  result  in  failure  or  permanent  set  in  the  propellant, 
generally  occurs  at  temperatures  higher  than  cure  temperature  after  sus¬ 
tained  exposure.  The  end  result  of  the  high  temperature  aging  process  may 
be  either  a  completely  degraded  propellant  or  an  essentially  new  propellant 
characterized  by  a  new  relaxation  spectrum,  higher  modulus,  and  reduced 
low  temperature  elongation. 
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3.4  CONSTANT  STRAIN  PROPERTIES 


As  indicated  in  Section  1.4,  the  constant  strain  test  as  commonly  employed  is 
diagnostic  of  propellants  with  rapidly  decaying  resistance  to  failure  under 
sustained  extension.  In  such  instances,  it  is  necessary  to  define  the  asymtotic 
level  of  strain  that  can  be  sustained  indefinitely  without  failure.  When  this 
level  is  more  than  2  percent  strain  below  the  constant  rate  strain  at  nominal 
maximum  stress  «m,  it  supercedes  «m  as  a  failure  criterion  for  the  evalu¬ 
ation  of  low  temperature  storage  and  firing  limits  of  a  motor. 

When  ballistic  necessity  dictates  the  application  of  a  propellant  inherently 
poor  in  constant  strain  properties,  it  is  necessary  to  evaluate  the  constant 
strain  asymtote  for  the  effects  of  temperature  and  initial  strain  rate  in  order 
to  fully  define  the  data  needed  in  the  motor  strain  analysis.  Problems  in  this 
category  require  specialized  treatment;  identification  of  the  existence  of  the 
problem,  however,  is  within  the  scope  of  standard  engineering  inquiry. 


3.5  PROPELLANT -LINER-SUBSTRATE  BOND  PROPERTIES 


The  bond-in-tension  apparatus  described  in  Section  1.5  provides  the  vehicle 
for  only  one  of  a  large  variety  of  bond  evaluation  methods  (see  the  ICRPG 
Solid  Propellant  Mechanical  Behavior  Manual,  1963,  for  description  of  ad¬ 
ditional  methods).  The  bond-in-tension  test  is  preferred  as  a  screening  test 
in  ordinary  development  and  quality  control  work  because  of  its  economy  in 
execution  and  the  ease  of  data  interpretation  on  a  limit  basis,  as  discussed 
in  Section  1.5.  In  general,  the  limit  analysis  simply  separates  excellent 
bonding  properties  (cohesive  propellant  failures)from  eminently  poor  bonding 
properties  (interfacial  failures)  and  isolates  inherently  weak  insulators  and 
liners  (cohesive  failure  in  liner  or  insulation).  Where  bond  strength  problems 
exist,  as  previously  noted,  a  chemical  or  process  problem  requiring  solution 
in  the  laboratory  is  indicated.  Where  optimization  of  case-bonding  strengths 
to  those  nearly  equal  to  the  observed  propellant  strength  is  not  possible,  ex¬ 
tremely  difficult  failure  analysis  problems  arise;  further  discussion  in  this 
connection  is  presented  in  Section  4. 
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3.6  REDUCED  LINEAR  COEFFICIENT  OF  EXPANSION 


The  propellant  coefficient  of  linear  thermal  expansion,  ®p,  is  ordinarily  a 
function  of  temperature  well  approximated  by 


«p  =  Ci  +  C2  T 


(3-42) 


where  Ci  and  C*  are  constants  and  T  is  temperature  (conveniently  T  in  °F  is 
used  to  define  Ci  as  the  value  of  <*p  at  0°F).  The _average  value  of  the  thermal 
coefficient  between  two  temperatures  Ti  and  T,  ffp,  is  given  by 


T 

»„  *  T7  -  TTTT  S  »p  iT.  (3-43) 

T  i 


where  is  the  specimen  length  at  Tj.  The  solution  for  thermal  strains  in  a 
case-bonded  motor  requires  integration  of  the  term  arp,  which  appears  in  the 
elastic  small  strain  solution,  or  equivalently,  direct  introduction  of  the  term 
«p  as  shown  in  Section  2.1.  The  simplest  procedure  for  generation  of  the 
thermal  coefficient  term  used  in  solution  for  thermal  strains  in  a  motor  is 
that  shown  in  Section  1.6.  The  technique  avoids  introduction  of  numerical 
error  through  the  redundant,  successive  differentiation  and  integration  of  the 
A-t/ L\  laboratory  measurement.  The  additional  simplification  introduced  by 
referencing  the  Af./ f-i  relationship  to  the  propellant  grain  zero  strain  temper¬ 
ature  permits  direct  correlation  of  laboratory  tensile  test  data  to  motor 
strains  in  failure  analysis. 
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3.7  BIAXIAL  TESTING  OF  TOXIC  OR  RARE  MATERIALS 


Rapid  evaluation  of  the  physical  properties  of  toxic  or  extremely  costly  pro¬ 
pellant  formulations  results  in  difficulties  in  executing  tensile  tests  because 
of  tensile  specimen  preparation  hazards  and/or  prohibitive  cost  requirements. 
The  diametral  compression  biaxial  test  has  been  used  as  a  means  for  amelio¬ 
rating  both  difficulties.  In  the  diametral  test,  a  circular  disc  is  compressed 
along  a  diameter,  producing  a  vertical  compression  and  a  lateral  tension,  at 
the  center  of  the  disc  (Figure  3-6).  In  extensive  tests,  it  has  been  found  that 
the  lateral  tensile  stress  and  strain  in  the  disc  agree  quite  accurately  to  the 
uniaxial  values  9m,  «m  (Section  1.2,  3.2).  The  diametral  test,  hence,  has 
been  used  as  an  auxiliary  method  for  screening  purposes  in  development  of 
toxic  and  expensive  formulations.  Specimens  consisting  of  pressed  discs 
J/4-inch  by  1 /4-inch  in  dimension  and  cast  discs  1  Vi-inch  by  Vi-inch  have  been 
used  successfully  in  this  method  of  evaluation.  In  both  types,  no  cutting  or 
trimming  of  the  specimen  was  required  and,  hence,  no  propellant  waste  in¬ 
curred.  In  addition,  because  the  physical  circumstances  of  the  test  facilitate 
the  collection  and  disposition  of  all  specimen  fragments,  cumulative  contami¬ 
nation  of  test  fixtures  with  toxic  materials  can  be  minimized.  The  following 
presents  the  test  procedure  used  in  conducting  and  analyzing  data  from  the 
tests. 


3.7.1  The  Diametral  Test 


Scope ;  The  following  describes  the  diametral  tension  -  compression  biaxial 
test,  presents  typical  test  data,  and  compares  it  to  uniaxial  tensile  data. 

Test  Apparatus:  Universal  tester  equipped  with  compression  platens  and 
load  cell,  and  with  temperature  conditioning  chamber. 

Test  Specimen;  A  disc  of  propellant  with  a  diameter -to -thickness  ratio  of 
approximately  3  is  preferred.  Diameters  of  6  inches  to  V-t-inch  have  been 
successfully  used.  Specimens  can  be  prepared  by  machining,  pressing  or 
casting,  but  regardless  of  preparation  method,  should  be  circular  within  1 
percent  and  opposing  planar  surfaces  should  be  parallel  within  1  percent. 

Test  Method  and  Analysis:  In  the  diametral  test,  a  right  circular  cylinder 
with  a  diame te r -to -thickne s s  ratio  of  approximately  3  of  the  material  to  be 
tested  is  placed  between  the  platens  of  a  compression  tester  (Figure  3-6). 
As  the  platens  move  together  along  a  diameter  of  the  specimen,  the  y-axis, 
a  resisting  force,  Fy,  is  developed  as  well  as  biaxial  stresses  9X  and  Oy  at 
the  center  of  the  sample. 

The  stress  in  line  with  the  compressive  force,  Fv,  is  o„  and  given  by  (Fitz¬ 
gerald,  1960),  y  y 
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Figure  3-6  Diametral  Compression  Test 


d  s  sample  diameter 
t  =  thickness 

the  negative  sign  signifies  a  compressive  stress. 

Along  the  axis  normal  tothe  applied  force  there  is  developed  a  tensile  stress, 
a  ,  given  by 


a 

x 


2F 

TTt  <4 


(3-45) 


the  corresponding  horizontal  strain  at  the  center  of  the  diametral  specimen 
is  given  by 


2Fy  (1  +  3v) 
*x  =  mfiEE 


(3-46) 


where  v  is  Poisson's  ratio  and  E  the  modulus  of  the  material. 
Vertical  strain  at  the  center  of  the  samples  is  given  by 


2F  +  (3  +  v) 
tt'cR'E - 


(3-47) 


The  total  deformation  along  the  center  horizontal  axis  U  (x,  o)  is  given  by 


u  (x*  o)  =  (  (1  -  v)  (2  -  tt)  +  2(1  +  v)  ] 


(3-48) 


The  ratio  of  horizontal  strain  at  the  center  of  the  sample  to  the  total  hori¬ 
zontal  unit  deformation  is  then 


x  =  2  (1  4  3  v) 

TJ/d  (1  -  v)  (2  -  tt)  +  (1  +  v) 


(3-49) 


This  ratio  of 


TT/d 


is  relatively  insensitive  to  Poisson's  ratio,  varying  from 


2.058  to  2.11  for  values  of  Poisson's  ratio  of  0.-3  and  0.5,  respectively. 
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The  test  is  carried  out  at  a  constant  crosshead  speed.  Measurements  of  the 
force  Fy  and  total  diameter  change  TT  are  made  (Figure  3-6  illustrates  one 
convenient  means  for  measuring  U)  and  the  stress  ox  and  strain  «x  at  the 
center  of  the  specimen  are  calculated  for  the  time  at  which  failure  by  cracking 
near  the  center  of  the  specimen  is  visually  observed. 

Biaxial  Test  Data  and  Discussion:  Data  for  diametral  compression  tests  of 
PBAA  propellant  conducted  at  -50®F,  100°F,  and  130°F  are  shown  in  Figure 
3-7.  Failures  at  0°F  and  lower  were  characterized  by  sharp  vertical  cracks 
appearing  along  the  vertical  axis  of  the  specimens.  Samples  at  30°F  and  70°F 
failed  along  lines  at  45  degrees  to  the  vertical  axis  (lines  of  maximum  con¬ 
stant  shear)  and  diametral  rupture  at  100°F  and  1 30°F  occurred  along  the 
vertical  axis. 

Strains  were  measured  photographically  at  the  center  of  the  specimens  and 
also  calculated  using  Equation  (3-49)and  the  total  horizontal  diameter  growth. 
Strains  calculated  thus  from  total  diametral  deflection  were  identical  within 
experimental  error  to  the  measured  strains  for  all  values  up  to  rupture. 

Failure  stress  values  calculated  for  the  center  of  the  specimen  using 

°x  (0'  y)  =  173?  (3"50) 


are  shown  in  Figure  3-7  with  a  failure-stress  curve  obtained  from  uniaxial 
tests.  It  should  be  noted  that  excellent  agreement  between  biaxial  failure 
stresses  and  uniaxial  maximum  stresses  was  obtained  even  through  biaxial 
stresses  were  calculated  from  infinitesimal  strain  theory.  A  few  percent 
deviation  only  was  found  between  uniaxial  strain  at  maximum  stress  and 
biaxial  strain  at  failure  values. 
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3.8  ANALOG  MOTOR  TESTS1 


Analog  motors  (or  strain  evaluation  cylinder s)  are  circular  port,  case-bonded 
motors  scaled  to  duplicate  the  strain-ver aus-temperature  profile  experienced 
in  a  flight  motor  of  circular  or  star -port  design.  Analog  motor  sizes  ranging 
from  2-inch  diameter  to  3-foot  diameter  have  been  used  in  simulation  of 
strains  in  flight  motors  ranging  from  5 -inch  to  1 20-inch  diameter s.  Four- 
and  seven-inch  diameter  analog  sizes  (Figure  3-8)  have  been  used  extensively 
at  LPC  for  the  following  reasons.  Four-inch  diameter  motors  are  the  small¬ 
est  size  in  which  sufficiently  accurate  measurements  of  grain  dimensions  can 
be  made.  While  most  of  the  4-  and  7 -inch  analogs  are  cast  in  nonfireable 
sleeves,  casting  and  ballistic  test  hardware  are  available  for  both  sizes  of 
motors.  The  availability  of  casting  tooling  decreases  the  cost  of  analog 
fabrication  and  the  availability  of  firing  tooling  facilitates  the  simultaneous 
investigation  of  ballistic  and  structural  behavior  under  controlled  conditions. 

Design  of  scaled  circular  port  grains  to  simulate  star-perforated  motors 
proceeds  from  consideration  of  the  thermal  strain  solution  presented  in  Sec¬ 
tion  2.  1.  The  pertinent  equations  are: 


where 


e 


£T 


*3 

(3-51) 

1  +  4h(1+Xj) 

[(ap  -  t'  T  a C)X2  -  Op]  P 

(3-52) 

%  K‘ 

(3-53) 

hoop  strain  at  port,  circular  port 

reduced  thermal  coefficient  of  grain  propellant 

temperature  increment  below  zero  strain  temperature 
(Section  2.  1) 


vc  =  case  Poisson's  ratio 

ac  =  case  thermal  expansion  coefficient 

X  =  b/ a  =  grain  O.  D.  /l.  D.  ratio 

P  =  Parr  end  condition  correction  for  finite  length  grains 
(Section  2.  1) 

=  star  valley  stress  concentration  factor  (Section  2.  1) 

‘Analog  test  techniques  are  also  discussed  in  detail  in  LPC-578-F-l,and  in 
the  ICRPG  Solid  Propellant  Mechanical  Behavior  Manual,  Section  4.  8.  Also 
see:  Ignatowski,  19o3. 
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Figure  3-8  Four  and  Seven -Inch  Diameter  Analog  Motors 
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A  circular  port  grain  scales  linearly  to  analog  sizes  by  keeping  equal  values 
of  X  and  length  to  diameter  (equal  F).  In  scaling  of  star  grains  for  circular 
port  analog  strain  simulation,  the  values  of  F  and/or  X  for  the  circular  port 
analog  grain  are  increased  to  compensate  for  the  K:  term  for  the  star  grain. 
P  increases  with  grain  length  to  diameter  ratio  and  decreases  with  increasing 
X  ;  hence,  some  juggling  of  these  values  is  necessary  to  arrive  at  the  dimen¬ 
sional  requirements  for  the  analog  design.  Generally,  it  is  convenient  to  take 
the  available  lengths  of  analog  motor  and  scale  X  appropriately  to  obtain  the 
required  strain  simulation. 

Parenthetically,  it  should  be  noted  that  analog  scaling,  keeping  F  for  the 
analog  equal  to  F  of  the  simulated  motor,  is  rigorously  more  desirable  since 
the  biaxial  stresses  in  the  analog  will  more  closely  simulate  those  of  the 
simulated  motor.  With  the  exception  of  motors  having  F  >0.99,  the  stress 
biaxiality  does  not  appear  to  be  particularly  significant.  When  F  for  a  motor 
is  equal  to  or  in  excess  of  0.99,  scaling  of  the  analog  simulation  motor  to 
maintain  P  of  the  analog  equal  to  or  in  excess  of  0.99  is  advisable. 

Analog  motors  are  conventionally  tested  by  controlled  thermal  cycling,  intro¬ 
ducing  heat  flow  only  through  the  case  wall  and  maintaining  the  grain  under 
desiccated  conditions.  A  typical  experimental  arrangement  is  shown  in  Fig¬ 
ure  3-9.  When  simulation  of  behavior  of  the  large  motor  conditions  is  desired, 
the  analog  cooling  (heating)  rates  are  controlled  such  that  the  strain-time 
variation  at  the  port  of  the  analog  reproduces  that  observed  for  the  large 
motor.  The  resulting  low  temperature  cracking  observations  in  the  analog 
will  generally  prove  to  be  slightly  conservative  (on  the  high  side)  with  respect 
to  the  full-scale  motor. 

When  analog  motors  are  used  as  vehicles  for  supplemental  laboratory  char¬ 
acterization  of  propellants  per  se,  shock  cycling,  extended  storage  at  low 
temperature,  etc.  ,  may  be  used  and  the  data  correlated  through  strain  rate 
and  temperature  correlation  where  the  rate  and  temperature  of  the  port  (at 
the  midpoint  of  the  grain)  are  the  parameters  of  interest. 

In  tests  of  the  analog  motors,  the  usual  observations  are  the  grain  port  di¬ 
ameter  at  the  grain  midpoint  (Figure  3-10)  visual  inspection  of  the  grain,  and 
temperature  measurements.  The  strain  at  the  grain  midpoint  is  that  pre¬ 
dicted  by  Equations  3-51  and  3-52  and  is  given  experimentally  by 


where  di  is  the  port  diameter  at  the  zero  strain  temperature  (or  the  mandrel 
O.D.)  and  Ad  is  the  change  in  port  diameter  with  respect  to  that  temperature. 
Visualinspection  of  analog  grains  with  the  unassisted  eye  is  inadvisable  since 
the  incipient  fissuring  of  grains  that  usually  precedes  crack  propagation  is 
quite  difficult  to  see.  Borescope  observation  is  usually  a  necessity  (Figure 
3-11).  Thermocouple  grain  and  case  temperature  measurement  then  com¬ 
pletes  the  required  set  of  experimental  data  necessary  for  correlation  of  rate 
and  temperature  behavior  of  the  analogs. 


3-32 


LOCKHEED  PROPULSION  COMPANY 


3-33 


9061 


LOCKHEED  PROPULSION  COMPANY 


Figure  3-9  Analog  Motor  Experimental  Set-Up 
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Figure  3-11  Borescope  Inspection  of  Analog  Motor  Grains 


It  should  be  noted  that  the  zero  strain  temperature  T,  appear s  to  be  influenced 
somewhat  by  analog  size,  particularly  for  analogs  smaller  than 4  by  14  inches 
cast  with  rubber-base  propellants.  A  method  for  independent  determination 
of  T,  is  duscussed  in  LPC  556-F - 1. 

In  summary,  analog  testing  can  be  regarded  as  second  only  to  full-scale 
evaluations.  Examples  of  the  results  of  analog  motor  tests  are  discussed  in 
LPC  578-F-l  of  this  report  and  by  Jones,  1963. 


1 
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Section  4 


GRAIN  STRUCTURAL  ANALYSIS  METHODS  DISCUSSION, 
ALTERNATE  METHODS,  AND  LIMIT  DESIGN  CONSIDERATIONS 


Background  discussion  and  additional  information  pertinent  to  the  engineering 
methods  described  in  Section  2  are  presented.  Limit  design  analysis  tech¬ 
niques  for  storage  and  launch  grain  stresses,  grain-case  terminations  and 
horizontal  star  point  deflections  are  described  in  Sections  4.6,  4.7,  and  4.8. 


4.  1  THERMAL  STRAIN  ANALYSIS 


The  thermal  strain  analysis  solution  given  in  Section  2.  I  considers  a  case- 
bonded  circular  port  grain  under  plane  strain  end  constraint  for  a  condition 
of  thermal  equilibrium.  Using  data  developed  by  Parr  (1962),  the  influence 
of  finite  grain  length  is  accounted  for.  These  data  are  reproduced  in  graph¬ 
ical  form  in  Appendix  D  with  the  permission  of  the  authors.  Using  the  data  of 
Fourney  and  Parmcrter  (1961),  the  perturbation  of  the  strain  field  by  geo¬ 
metrically  regular  star  port  perforations  is  accounted  for.  These  data  also 
are  reproduced  in  Appendix  D  with  the  permission  of  the  authors. 

The  equivalent  iteration  technique,  which  compensates  for  finite  grain  port 
deflection  (Reference  Note  2.  11,  Section  2.  l),  is  an  empirically  developed 
operation  that  has  been  found  to  be  within  one  per  cent  agreement  with  precise 
computer  iteration  solutions  at  strain  levels  up  to  25  percent.  The  neglect  of 
thermally  induced  shear  stresses  near  the  grain-case  boundary  produces  a 
trivial  effect  upon  the  predicted  strains.  The  detailed  analysis  for  the  plane 
strain  cylinder  is  included  in  Williams,  et.  al,  1961  (GALCIT  SM  61-5);  the 
treatments  for  end  conditions  and  internal  star  geometry  are  presented  in  the 
previously  referenced  publications  of  Parr,  Fourney  and  Parmerter.  Addi¬ 
tional  details  of  the  thermal  strain  analysis  problem  and  experimental  data 
are  included  in  LPC  578-F-l  and  LPC  556-F-l.  Further  discussion  herein 
is  limited  to  the  key  factors  that  require  emphasis  in  engineering  applications. 

The  given  solution  predicts  the  tensile  hoop  strain  at  the  midpoint  of  a  grain 
for  a  propellant  of  known  coefficient  of  thermal  expansion  and  zero  strain 
temperature  T,  (Section  1.6).  Grain  hoop  strain  is  ordinarily  a  maximum  at 
this  location  and  localizes  web  crack  initiation  there.  Deviations  from  ther¬ 
mal  equilibrium  in  the  grain  produce  errors  in  strain  prediction:  upper  and 
lower  limits  for  the  strains  can  be  calculated  in  such  instances  by  taking 
equilibrium  strain-ver sus -temperature  data  for  the  upper  and  lower  bounds 
of  the  grain  temperature.  The  strain  at  thermal  equilibrium  at  a  given  tem¬ 
perature  will  be  the  upper  limit  of  strain  at  that  temperature. 
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Referencing  the  strain  calculations  to  the  temperature  T,  is  necessary  in  pre¬ 
dicting  grain  failure,  as  discussed  in  Section  4.  2.  The  typical  values  given 
for  grain  zero  strain  temperatures  lead  to  reasonably  accurate  grain  failure 
analysis.  An  experimental  technique  for  independent  measurement  of  T l  is 
given  in  LPC  556-F-l.  In  the  absence  of  other  experimental  information, 
determination  of  the  temperature  at  which  a  grain  I.  D.  is  equal  to  the  cor¬ 
responding  mandrel  O.  D.  yields  a  reasonable  estimate  of  the  zero  strain 
temperature.  It  may  be  observed  that  long  term  storage  or  cycling  of  a 
motor  produces  changes  in  the  grain  dimensions  in  addition  to  changes  caused 
purely  by  thermal  contraction  effects.  In  general,  changes  in  dimension 
that  tend  to  reduce  the  grain  port  I.  D.  (as  measured  at  a  given  temperature) 
are  indicative  of  incipient  failure  in  the  grain,  arising  possibly  from  binder- 
filler  dewetting.  Increases  in  port  I.  D.  are  indicative  of  additional  pro¬ 
pellant  cure  or,  at  low  temperature,  binder  crystallization.  In  the  latter 
instances,  the  properties  of  the  propellant  are  also  subject  to  severe  alter¬ 
ation  and  special  methods  are  required  in  structural  analysis  of  the  grain. 
The  zero  strain  temperature  of  the  grain  can  be  expected  to  change  in  these 
cases,  invalidating  the  grain  failure  calculations  to  an  unknown  extent.  In 
fully  cured  PBAA,  Polysulfide  and  PBAN  formulations  such  changes  are 
usually  of  minor  extent. 

The  analysis  as  given  assumes  a  propellant  Poisson's  ratio  of  ft,  i.  e. ,  the 
usual  assumption  of  rubber- like  incompressibility.  While  a  certain  amount  of 
controversy  exists  with  regard  to  the  assumption,  the  following  summarizes 
pragmatic  observations  in  this  regard.  Thermally  induced  grain  strains  have 
been  accurately  predicted  in  motor  sizes  ranging  from  four-inch  to  three- foot 
diameters  at  temperatures  between  110° and  +200°F  in  more  than  one  hundred 
motors  examined  under  conditions  of  rigorously  controlled  environment  (LPC 
556-F-l;  LPC  578-F-l;  Jones,  1963).  The  Poisson's  ratio  for  the  propel¬ 
lants,  which  ranged  up  to  86  percent  solids  loading  by  weight,  was  taken  as 
7j>.  Measured  port  strain  levels  ranged  up  to  25  percent.  In  all  such  cases 
examined,  grain  measurements  indicating  apparent  values  of  Poisson's  ratio 
of  less  than  '/i  exhibited  grain  failure,  in  the  majority  of  cases  by  web  crack¬ 
ing.  Strains  in  sound  grains  were  predicted  analytically  to  within  a  percent 
of  strain.  Hence,  the  assumption  is  considered  to  be  accurate  for  the  con¬ 
dition  of  interest  and  cannot,  by  itself,  lead  to  other  than  conservatively  high 
strain  predictions. 

The  Parr  end  condition  correction  F  was  developed  using  infinitesimal  elastic 
theory  and  numerical  data  developed  by  computer  calculation  at  Rohm  and 
Haas  Redstone.  Within  the  assumption  of  infinitesimal  deformations,  previ¬ 
ously  unpublished  experimental  data  at  LPC  have  shown  the  values  of  P  to  be 
exact  within  the  limits  of  experimental  accuracy. 

The  star  valley  stress  (strain)  concentration  factor  Kj  data  as  developed  by 
Fourney  and  Parmerter  also  considered  only  infinitesimal  strain  conditions. 
Generally,  thermal  cooling  produces  an  increase  in  the  radius  of  star  valley 
curvature,  tending  to  reduce  the  stress  concentration  factor.  Values  of 
based  on  zero  strain  temperature  grain  (or  mandrel)  dimensions  lead  in  any 
case  to  slight  conservatism  in  that  predicted  strains  will  tend  to  be  high.  In 
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instances  where  the  star  valley  geometries  deviate  from  a  circular  arc.  esti¬ 
mates  of  Ki  using  the  data  of  Fourney  and  Parmerter  are  inaccurate.  While 
bounding  values  may  be  estimated  using  minimum  and  maximum  radii  for  the 
star  valley  geometry,  it  is  usually  necessary  to  obtain  values  for  Kj  by  photo¬ 
elastic  modeling  of  the  actual  grain  cross-section.  The  techniques  required 
in  experimental  measurement  of  are  specialized  in  nature;  Fourney  and 
Parmerter,  1961,  present  an  up  to  date  discussion  of  the  techniques  and 
pertinent  source  references. 


4.2  THERMAL  STRAIN  FAILURE  ANALYSIS 


The  method  for  grain  failure  prediction  presented  in  Section  2.  2  is  an  engi¬ 
neering  approximation  to  a  cumulative  strain  failure  criterion  proposed  by 
Williams,  et.  al,  1961  (GALCIT  SM  61-5)  and  discussed  by  Jones,  1963.  For 
the  particular  condition  of  interest,  monotonic  cooling  to  failure,  cumulative 
strain  failure  calculations  agree  with  the  method  shown.  In  spite  of  the  fact 
that  thermal  gradient  conditions  and  attendant  port  hoop  strain  rates  are 
treated  in  apparently  cavalier  fashion,  the  technique  works  quite  accurately, 
(LPC  578  -F-l;  Jones,  1963  ).  The  dominating  factor  in  determining  low 
temperature  grain  cracking  is  the  rapid  decrease  in  propellant  strain  capa¬ 
bility  at  low  temperatures.  At  temperatures  near  the  glass  transition  region 
(approximately  -120°F  in  the  example  shown  in  Section  2.  1  and  2.2),  the 
properties  of  the  propellant  vary  about  equally  for  a  ten-degree  change  in 
temperature  or  a  one  order-of-magnitude  change  in  strain  rate;  hence,  the 
relative  unconcern  with  precision  in  specifying  strain  rates  in  this  region. 
In  consideration  of  the  effects  of  thermal  gradients  it  should  be  noted  that  in 
cooling  to  a  given  temperature,  the  thermal  gradient  disappears  as  that 
temperature  is  approached.  Because  of  the  greater  sensitivity  of  the  propel¬ 
lant  properties  to  temperature  than  to  strain  rate  in  the  low  temperature  re¬ 
gime,  the  final  stages  of  cooling  are  of  greatest  significance  in  determining 
failure.  The  physical  conditions  in  a  grain  approaching  thermal  equilibrium 
near  its  cracking  limit  dictate  a  convergence  toward  the  assumed  conditions 
of  thermal  equilibrium  and,  thus,  a  convergence  in  the  accuracy  of  failure 
prediction.  Errors  introduced  from  this  cause  tend  to  lead  to  prediction  of 
conservatively  high  grain  cracking  temperatures.  In  short,  the  procedure 
leads  to  prediction  of  minimum  safety  factors  for  the  worst  condition  of  ther¬ 
mal  cooling  ordinarily  imposed  on  a  motor. 

Where  more  complicated  conditions  of  cooling  require  analysis,  the  procedure 
is  as  follows.  Starting  from  the  zero-time  zero  strain  condition,  the  temper¬ 
ature,  strain  and  strain  rate  at  the  port  of  a  grain  are  calculated  incre¬ 
mentally.  In  a  small  temperature  increment,  the  ratio  of  the  incremental 
strain  to  the  strain  capability  in  the  temperature  increment  and  at  the  mean 
strain  rate  is  calculated  and  taken  to  be  the  fractional  use  of  the  strain  capa¬ 
bility  of  the  propellant.  For  conditions  of  monotonic  cooling,  then,  failure  iB 
predicted  when  the  sum  of  the  fractional  usage  of  strain  with  an  increasing 
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number  of  temperature  increments  reaches  the  value  unity.  This  technique 
requires  considerable  precision  in  all  phases  of  the  analysis  and  is  presently 
regarded  to  be  a  matter  of  concern  for  specialists  in  grain  analysis. 

Comment  is  required  on  the  correlation  of  propellant  uniaxial  failure  proper¬ 
ties  to  prediction  of  motor  failure  under  the  multiaxial  strain  field  existing 
at  the  port  of  either  circularly  or  star  perforated  grains.  Laboratory  multi¬ 
axial  failure  data  suggests  that  propellant  strain  capabilities  in  strong  triaxial 
tensile  strain  fields  are  much  smaller  than  in  uniaxial  tension.  Since  the  ther¬ 
mal  strain  field  in  a  cooled  grain  tends  to  be  triaxially  tensile,  it  may  be 
suspected  that  uniaxial  data  would  predict  erroneously  low  motor -cracking 
temperature.  Noting  that  intensive  research  in  this  connection  is  currently 
in  progress  throughout  the  propellant  industry,  the  present  status  of  engi¬ 
neering  knowledge  in  this  connection  is  best  summarized  as  follows. 

Properly  and  precisely  executed,  the  simplified  engineering  analyses  pre¬ 
sented  in  this  handbook  have  been  demonstrated  to  be  either  accurate  within 
the  accuracy  limits  stated  or  conservative  in  all  instances  where  either  model 
or  flight  motors  have  been  cooled  to  failure  at  LPC  during  the  past  two  years. 
The  experience  cited  includes  propellants  basedon  PBAN,  Polysulfide,  PBAA, 
and  terminally  carboxylated  binders  in  motors  ranging  in  size  from  4-inch  to 
120-inch  diameters.  Proper  analysis  of  the  laboratory  physical  property 
data  is,  of  course,  vital  to  the  analysis  (Reference  Section  1  and  3). 

Referencing  the  grain  strain  calculations  to  the  zero  strain  temperature  is 
essential  because  the  total  strain  imposed  on  the  propellant  regulates  the 
initiation  of  failure  and  is  indeed  the  strain  measured  in  the  laboratory  tensile 
test.  Omission  of  the  calculation  step  in  which  the  subtractive  term  a  is 
neglected  (Reference  Table  2-1  and  Note  2.16),  is  also  essential  in  correlating 
failure  strains  in  the  grain  to  the  corresponding  strains  as  measured  in  the 
laboratory. 
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4.3  GRAIN  SLUMP  DEFLECTIONS 


The  solutions  for  slump  deflections  in  circular  port  grains  subjected  to  axial 
accelerations  given  in  Section  2.3  were  extracted  from  Parr,  1963,  with  the 
author's  permission.  The  numerical  methods  used  by  Parr  in  obtaining  the 
dimensionless  form  of  the  pertinent  solutions  are  discussed  in  the  above 
referenced  report  and  in  Parr,  1962. 

In  each  of  the  cases  considered  by  Parr,  an  elastic  solution  was  obtained  re¬ 
lating  a  deflection  U  to  the  imposed  acceleration  ng,  the  density  of  the  propel¬ 
lant  P,  the  radial  dimension  of  the  case  b,  the  propellant  (elastic)  modulus  E, 
and  a  numerical  constant  p  which  is  dependent  on  the  grain  web  fraction,  con¬ 
straint,  and  length  to  outside  diameter  ratio: 

U  =  P  —  P^S. 


The  motor  case  was  considered  to  be  rigid  in  all  instances.  The  viscoelastic 
solution  for  the  time  dependent  deflection  U(t)  is  obtained,  in  terms  of  the 
differential  operator  notation  (Section  3.1)  as 

UCt)  -  §  (Pb2  p  J  n-g(t) 

where  P  and  Q  are  linear  differential  operators  and  the  notation  n*g(t)  indi¬ 
cates  the  time  dependence  of  the  applied  axial  acceleration.  Rigorously,  the 
solution  for  a  step  function  acceleration 

t  n-g(t)  =  0,  t  <_  0 
n-g(t)  =  n- g  U(t),  ] 

‘  n-g(t)  =  n.  g,  t  >  0 

is 


*  n  -(3.  t 

U(t)  =  Pb2  P  n- g  .2  B.  e  1  , 

1  =  0  1 

where  the  notation  (Section  3.1  and  Appendix  B)  indicates  the  compliance  form 
of  the  propellant  viscoelastic  stress-strain  law.  In  approximating  this  so¬ 
lution  for  engineering  calculations  the  nearlyexact  relationship  betweencreep 
compliance  D(t)  and  relaxation  modulus  E(t),  in  uniaxial  tension, 

E(t)  =  W) 
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was  introduced.  Hence, 


U(t)  =  £b*  P  n-g  D(t) 


*  * 

~  Pbz  p  n* 

=  — ETtr 


the  latter  form  of  the  equation  being  the  one  on  which  the  numerical  solutions 
shown  in  Section  2.3  were  based.  The  solutions  thus  are  approximations  that 
are  estimated  to  be  within  tolerable  engineering  limits  of  accuracy,  with  re¬ 
spect  to  the  viscoelastic  assumptions.  The  accuracy  of  the  elastic  solutions 
of  Parr,  insofar  as  experimental  data  are  available,  is  also  within  the  limits 
of  engineering  accuracy  (LPC  556-F-l,  Fitzgerald,  et.  al,  1961). 


4-7 


LOCKHEED  PROPULSION  COMPANY 


4.  4  GRAIN  PRESSURIZATION 


The  limit  solutions  for  grain  pressurization  strains  were  obtained  from  the 
Lame'  solutions  for  stresses  in  a  right  circular  cylinder  (for  example,  see 
pp.  125-145,  GALCIT  SM  61-5).  The  limit  solution  for  end-bonded  grains 
is  the  plane  strain  analysis  for  a  case-bonded  grain  as  given  in  the  cited 
reference  for  the  condition  where  the  modulus  of  the  propellant  grain  is  taken 
to  be  zero.  The  solution  considers  a  radially  expanding,  longitudinally  rigid 
motor  case.  In  the  limit  solution  given  in  Section  2.  4,  the  propellant  is 
assumed  to  be  incompressible  (i.  e.  ,  Poisson's  ratio  is  taken  as  1/ 2)  and  for 
the  case  Poisson's  ratio  is  taken  as  0.3.  The  limit  solution  for  the  grain  with 
unbonded  ends  proceeds  with  the  assumptions  as  given  above  except  that  the 
grain  longitudinal  stress  c2  is  specified  as  oz  =  -pj,  where  p^  is  the  motor 
pressure.  In  both  solutions  only  the  strain  at  the  port,  the  largest  strain,  is 
given. 


The  grain  strains  at  the  port  in  the  cases  considered  increase  from  zero 
toward  the  limits  given  in  Section  2.  4  with  a  strain-time  variation  controlled 
by  the  viscoelastic  properties  of  the  grain  and  the  time  history  of  pressuriza¬ 
tion.  Viscoelastic  solutions  for  the  strains  are  not  usually  pertinent  for  mod¬ 
erate  web  fraction  grains  in  metal  cases  since  the  limit  solutions  show  trivial 
limiting  values  for  the  strains.  Instances  where  the  pressurization  limit 
strains  exceed  two  or  three  percent  present  special  analysis  and  design  prob¬ 
lems  beyond  the  present  scope  of  this  handbook. 


The  effect  of  finite  compressibility  on  the  pressurization  strains,  by  limit 
analysis,  is  given  by 


€ 


0  p  ' 


(X2-l) 

2 


Pi 

TC 


for  the  plane  strain  (ends  bonded)  condition  where  K  is  the  propellant  bulk 
modulus,  p:  the  internal  pressure  and  \  the  O.  D.  to  I.  D.  ratio  for  the  grain. 
Typical  values  for  the  contribution  of  finite  grain  compressibility  to  grain 
pressurization  are  0.3  and  1.5  percent  for  50  and  75  percent  web  fraction 
grains  operating  at  1000  psi  internal  pressure  for  a  typical  value  of  K  =  5x  105 
psi/in.J.  Grainhoop  strains  from  compressibility  effects  are  lower  for  grains 
with  unbonded  ends. 


In  summary,  grain  pressurization  strains  in  moderately  designed,  steel-case 
rocket  motors  are  generally  of  minor  importance,  except  very  near  the  low 
temperature  storage  cracking  limit  of  the  grain.  In  this  region  and  in  the 
instances  of  high  web  fraction  high  pressure  operation  and/or  low  modulus 
casings,  highly  specialized  analytical  techniques  are  required  in  precise 
evaluation  of  grain  capabilities. 
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4.  5  PRESSURIZATION  FAILURE  ANALYSIS 


The  cumulative  damage  criterion  applied  in  pressurization  failure  analysis 
has  not  as  yet  been  adequately  demonstrated  for  solid  propellants.  Hence,  its 
adoption  for  this  purpose  is  pragmatic  in  nature  in  that  it  is  consistent  with 
engineering  practice  and  has  not  been  contradicted  by  experimental  motor 
firings  at  LPC.  Research  in  cumulative  damage  processes  is  presently  in 
progress  at  LPC. 

Selection  of  the  port  hoop  strain  as  the  failure -producing  phenomenon  of 
direct  interest  is  founded  on  both  theoretical  and  experimental  grounds. 
Viscoelastic  analyses  for  hoop  stresses  in  dynamically  pressurized  case- 
bonded  propellant  grains  show  that  few  if  any  modern  propellants  possess 
transient  modulus  values  sufficiently  large  to  permit  the  development  of  hoop 
tensile  stress  at  the  port  except,  possibly,  when  the  grain  temperature  is  so 
low  that  glassy  behavior  may  be  expected.  In  general,  grain  failure  from 
thermally  induced  strains  will  be  predicted  at  higher  temperatures  in  such 
instances.  The  hoop  stresses  in  a  grain  can  be  expected  to  be  compressive 
in  most  instances  and  therefore  cannot  produce  stress- determined  tensile 
failure.  The  failure  of. propellants  in  a  state  of  tensile  strain  and  hydrostatic 
compression  has  been  shown  to  be  relatively  independent  of  the  superposed 
hydrostatic  field,  and  in  reasonable  agreement  with  the  tensile  failure  strains 
as  measured  by  standard  techniques  (Kruse  and  Mahaffey,  1963). 
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4.6  STORAGE  AND  LAUNCH  STRESSES 


Shear  stresses  at  the  propellant-case  interfacial  area  in  vertically  stored  or 
axially  accelerated  grains  have  been  a  matter  of  concern,  particularly  in  the 
vicinity  of  the  top  grain  termination  circumference.  In  general,  the  geometric 
conditions  at  the  termination  circumference  preclude  a  precise  analytic  evalu¬ 
ation  of  stress  conditions  there,  and  it  has  been  expedient  to  proceed  in  prag¬ 
matic  consideration  of  the  problem.  Present  design  practice  at  LPC  requires 
that  the  average  shear  stress  xz(b)  at  the  case  be  <76  the  tensile  strength  of 
the  propellant,  as  compared  in  the  appropriate  viscoelastic  time  and  tem¬ 
perature  envelope.  The  stress  Tz(b)  is  defined  as 


Tz(b)  = 


where 


n -g  =  acceleration  in  gravities 
W  =  total  weight  of  the  grain 
A  =  bonded  area  of  the  grain 

With  reservations,  the  analysis  is  applied  to  grains  independently  of  the 
geometry  of  the  grain  end,  so  long  as  a  stress  relief  boot  or  flap  terminates 
the  grain  at  the  motor  case  wall  (Ref.  Section  4-7).  Only  unpressurized 
motors  are  considered:  in  pressurized  motors  the  bond  stresses  are  sub¬ 
jected  to  a  superposed  compressive  stress  that  compensates  for  the  effects  of 
the  acceleration  body  forces  except,  possibly,  in  extremely  high-acceleration 
motors. 

The  rationale  for  the  design  relationship  assumes  a  1  to  2  relationship  between 
a  shear  stress  and  its  uniaxial  tensile  equivalent  and  accounts  for  a  stress 
concentration  at  the  termination  of  3.  The  relationship  is  considered  to  be 
conservative,  and  the  fact  that  it  is  conservative  has  been  borne  out  in  a 
limited  number  of  experiments.  Research  in  this  connection  is  presently  in 
progress  at  LPC. 

In  illustration  of  the  Tz(b)  <p/ 6  relationship  in  design  evaluation,  consider  a 
50  percent  web  fraction,  240-inch  diameter  circular  port  case -bonded  grain 
with  unsupported  ends  stored  with  its  axis  vertical.  The  shear  stress  t  (b) 
is,  approximately,  3*4  psi  and  the  equivalent  tension  21  psi.  The  propellant- 
case  bond  constant  load  strength  (see  Section  1.  5)  required  for  conservative 
design  for  a  storage  temperature  of  70°F  would  be  21  psi  at  the  maximum 
anticipated  storage  time.  Similarly,  unpressurized  and  under  a  launch  ac¬ 
celeration  of  10  gravities  for  two  minutes,  the  required  minimum  constant 
load  bond  strength  would  be  210  psi  for  two  minutes. 

This  analysis  ispresentedas  a  reasonable  and  conservative  engineering  esti¬ 
mate,  not  as  a  necessarily  valid  means  for  specifying  grain  structural  safety 
factors.  1 
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4.7  GRAIN  TERMINATIONS 


One  of  the  most  common  grain  failures  is  the  circumferential  cracking  that 
occurs  at  the  grain-case  termination  location  (Figure  4-  1  ).  Such  cracks 
typically  propagate  along  the  propellant-case  interfacial  zone,  slightly  inside 
the  propellant  grain,  and  stop  permanently  at  some  distance  from  the  ends  of 
the  grain.  As  applied  to  specification  of  the  stresses  at  the  termination  cir¬ 
cumference,  elastic  analysis  predicts  infinite  stress  at  this  location  when  a 
geometric  discontinuity  exists.  The  prevalence  of  the  failures  experienced 
when  the  grain-case  te rmination  does  include  a  geometric  discontinuity  bears 
out  the  analysis  in  that  they  present  prima  facie  evidence  for  the  existence  of 
high  stresses  there. 

The  observed  behavior  of  the  failures,  i.  e.  ,  the  finite  propagation  depth,  has 
been  predicted  for  the  case  of  interest  in  an  unpublished  analysis  by  J.  E. 
Fitzgerald  based  on  a  published  analysis  of  M.  L.  Williams.  The  analysis 
predicts  that  the  depth  of  propagation  along  the  interface  should  vary  with 
motor  diameter  as  indicated  in  Figure  4-2,  which  also  shows  measured  crack 
depths.  The  specific  influence  of  other  motor  geometric  parameters,  such 
as  web  fraction,  irregular  case-end  geometry,  and  motor  length-to-diameter 
ratio  is  uncertain.  The  cited  information  applies  generally  to  motors  of 
approximately  50  percent  web  fraction,  length-to-diameter  ratio  1  to  3,  and 
of  cylindrical  design  as  indicated  in  Figure  4-1. 

In  the  relative  absence  of  precise  stress  and  deformation  field  design  infor¬ 
mation  relative  to  the  termination  circumference,  the  obvious  engineering 
approach  of  stress  relief  modification  of  the  geometry  has  been  found  to  be 
effective.  Typical  grain  terminations  modified  for  stress  relief  are  shown  in 
Figure  4-3  .  As  applicable  in  the  geometries,  the  stress  relief  flaps  extend 
to  approximately  the  depths  from  the  ends  of  the  grains  where  crack  propa¬ 
gation  would  be  expected  to  stop,  based  on  the  data  of  Figure  4-2  . 

The  selection  of  materials  for  the  flap  of  "boot"  is  governed  by  the  necessity 
for: 

(a)  Compliance  equal  to  or  greater  than  that  of  the  propellant  at  all 
anticipated  use  temperatures. 

(b)  Tensile  strength  and  elongation  greater  than  that  of  the  propel¬ 
lant  at  all  anticipated  use  temperatures. 

(c)  Insulation  capability  consistent  with  the  ballistic  environment. 

In  general,  most  high-quality,  silica-filled  insulation  rubbers  are  adequate 
for  the  purpose.  It  is  vital  that  thebondof  the  boot  to  its  usual  insulation  sub¬ 
strate  and  the  insulation  substrate -to-case  bond  both  be  strong  and  continuous . 

The  use  of  boot  relief  at  the  ends  of  motors,  as  shown  in  the  lower  examples 
of  Figure  4-3  ,  tends  to  minimize  thermal  strain  requirements  on  moderate 
length  grains  (Reference  Section  2.1,  4.1)  and  minimize  bonding  difficulties 
at  the  ends  of  grains  of  any  length. 
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Figure  4-1  Grain  Case  Termination  Cracking 
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Figure  4-2  Grain  Termination  Crack  Depth  versus  Motor  Diameter 
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Figure  4-3  Typical 


4.8  STAR  POINT  DEFLECTION  DURING  HORIZONTAL  STORAGE 


The  rigorous  analysis  for  star  point  slump  deflection  requires,  as  a  starting 
point,  a  solution  within  the  theory  of  elasticity  (Timoshenko  and  Goodier,  1951). 
However,  under  consideration  of  the  appropriate  field  equations  and  boundary 
conditions,  useful  information  can  be  obtained  by  employing  simple  beam 
analysis  for  an  estimate  (Timoshenko,  1956).  Although  a  star  point  is  not 
supported  in  a  simple  manner  near  the  motor  case  it  is  expedient  to  assume 
that  it  deflects  like  a  cantilever  beam  of  nonuniform  cross  section  which  is 
supported  rigidly  at  the  wall  (Figure  4-4  ).  The  star  point  maybe  approxi¬ 
mated  geometrically  either  by  a  truncated  or  nontruncated  triangular  section. 
If  the  star  point  cannot  be  represented  reasonably  by  such  a  configuration, 
the  following  equations  do  not  apply. 


4.8.1  Deflection  Expressions 


The  deflection  Y  at  a  star  point  tip  in  a  motor  in  horizontal  storage  (Figure 
4-  5  )  can  be  written  as 

Y  =  YB  +  V 

For  the  non-truncated  triangular  beam  (Figure  4-  5  )  one  obtains  for  the 
elastic  solution: 


v  _  3  P  L4 

YB  -  7  E  H7 

YS  =  K13  £  L* 

where: 

p  =  Propellant  density  (lb/ in.3) 

E  =  Propellant  modulus  (psi) 

L,  H  =  As  defined  in  Figure  4-  5  ,  in  inches. 

The  star  point  is  assumed  to  be  long  with  respect  to  the  dimension  L. 

If  the  star  point  is  represented  by  a  moderately  truncated  section  (Figure 
4-  6),  the  appropriate  expressions  are,  for  hQ/H  <0.1: 


+  1.13 
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Figure  4  -6  T run;at<?d  Star  P<.  irt  Geometry 
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1.13P  /  L  \* 

Ys  *  t  \irj 


S1  * 2  i?S  ■ 

where  p  and  £  are  as  defined  above  and  the  dimension  L,  H,  and  h  are  as 
shown  in  Figure  4-  6  . 

4.8.2  Time  Dependence  of  Slump  Deflection 


The  time  dependent  forms  of  the  preceding  equation  are  obtained  approxi¬ 
mately  by  replacement  of  E(t)  for  E,  where  the  numerical  values  of  E(t)  may 
be  obtained  for  the  appropriate  temperatures  from  Figure  1-  7  ,  Section  1.1. 
The  solutions  given  will  generally  tend  to  give  conservatively  large  estimates 
of  maximum  star  point  tip  deflections. 
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Section  A-l 


TEST  SPECIMEN  PREPARATION 


A -1.1  BACKGROUND 

The  essential  purpose  in  laboratory  testing  of  solid  propellants  is  to  deter¬ 
mine  the  physical  characteristics  of  the  propellants  as  they  exist  in  actual 
rocket  motors.  It  is  of  critical  importance,  therefore,  that  the  laboratory 
test  specimens  be  cut,  machined,  and  handled  using  techniques  that  do  not 
degrade  the  bulk  characteristics  of  the  propellant.  The  practical  problem  is 
one  of  imposing  the  minimum  clamping  and  handling  stresses  that  permit  the 
production  of  accurately  dimensioned  test  specimens.  A  second  problem,  but 
not  an  unimportant  one,  is  that  of  performing  the  necessary  specimen  prepa¬ 
ration  operations  while  maintaining  a  high  degree  of  personnel  protection 
against  the  possibility  of  fire  or  explosion.  The  following  technique  is  repre¬ 
sentative  of  the  recommended  procedures  for  preparation  of  JANAF  and 
bonded  tab  end  tensile  specimens. 


A-l. 2  PROPELLANT  CUTTING- -GENERAL 


Bulk  samples  of  propellant  for  reduction  to  te  st  specimen  form  are  customarily 
received  in  the  form  of  cylindrical  or  rectangular  blocks  ranging  in  volume 
size  between  one  quart  and  one  gallon.  Modified  bandsaws  are  used  for  cutting 
the  bulk  samples  to  various  section  shapes  slightly  larger  than  the  size  of  the 
desired  final  test  specimen;  final  sizing  and  shaping  of  specimens  are  done  by 
machining. 

A  typical  saw  (in  use  at  LPC)  is  an  industrial  two-wheel,  vertical  band  saw 
(Figure  A-l)  modified  for  use  in  propellant  cutting  as  described  below.  The 
conventional  electric  motor  was  replaced  with  an  explosion-proof  motor;  metal 
blade  guides  were  changed  to  teflon,  a  vacuum  scavenging  port  was  installed 
for  collection  of  small  cuttings  and  oxidizer  dust,  and  the  existing  sample 
vise  was  modified  for  gripping  the  various  sample  containers  present  in  pro¬ 
pellant  sample  cutting.  Conventional  skiptooth  bandsaw  blades  with  high-side, 
rake  angle  teeth  are  used  in  routine  cutting. 

In  cutting  propellant,  it  was  found  that  vacuum  scavenging  was  only  partially 
effective,  consequently  manual  cleanup  of  propellant  fragments  between  saw 
cuts  is  a  necessary  safety  precaution.  The  sawing  procedure  followed  is  one 
in  which  the  actual  saw  cut  is  made  using  remote  control  from  outside  the 
sawing  barricade  to  actuate  the  hydraulically  operated  table  feed  (Figure  A-2)^ 
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Figure  A-l  Vertical  Band  Saw  Adapted  for  Remote  Operation 
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Redaction  of  a  bulk  sample  takes  place  by  sawing  and  resawing  until  blocks, 
bars,  slabs,  etc.,  of  the  desired  rough  size  are  obtained.  The  modified  vise 
on  the  saw  permits  the  preparation  of  blocks  and  bars  with  cubic  or  rectangu¬ 
lar  dimensions  between  one-half  and  seven  inches.  The  vise  used  in  clamp¬ 
ing  samples  for  saw  cutting  is  a  conventional  screw  type  vise,  except  for 
elimination  of  metal-to-metal  contact. 

The  saw  installation  is  set  up  with  both  automatic  and  manual  water  drowning 
provisions.  Ordinarily,  the  saw  operator  can  operate  the  manual  valve  faster 
than  the  fusible  valves  at  the  saw  melt  when  a  propellant  fire  starts.  In  either 
case,  the  water  sprays  minimize  damage  to  the  equipment  when  fires  occur. 
The  vacuum  system  uses  only  plastic  hose  conduit  and  has  an  intermediate 
drowning  tank  for  collection  of  propellant  dust  at  the  saw. 


A  -  1 . 3  PROPELLAN  T  MILLING- -GENERAL 


Propellant  samples  rough  cut  by  sawing  are  machined  to  final  size  using  a 
modified  milling  machine  (Figure  A-3).  The  basic  modifications  to  the  com¬ 
mercial  model  machine  are  described  herein.  The  standard  electric  motors 
for  the  cutting  head  and  table  feed  drives  are  replaced  with  explosion-proof 
motors,  and  the  belt  drive  mechanism  in  the  head  of  the  machine  is  covered 
and  sealed.  A  sample  holding  fixture,  complete  with  dust  cover  and  cutting 
basket,  is  fixed  to  the  mill  table.  Vacuum  scavenging  arrangements  and  water 
drowning  heads  are  included  in  the  installation.  The  operational  features  of 
the  drowning  system  are  identical  to  those  of  the  saw  operation  previously 
described.  In  general,  the  methods  employed  with  the  mill  are  the  same  as 
those  used  in  saw  operation.  All  operations  that  must  be  done  manually  in 
the  mill  barricade  are  designed  so  that  metal-to-metal  frictional  contact  is 
eliminated  in  adjustment  of  the  mill  between  passes  and  clamping  of  specimens 
on  the  mill  table.  Manual  scavenging  of  propellant  cuttings  is  necessary  be - 
tween  cuts.  Operationof  the  mill  during  the  actual  cutting  ope  ration  is  carried 
out  by  remote  control. 

The  cutter  heads  used  in  all  milling  operations  are  of  the  circular  fly  cutter 
type  (Figure  A -4).  These  cutters  are  adjusted  variously  for  cutting  the  re¬ 
entrant  section  on  JANAF  tensile  specimens  or  for  milling  planar  surfaces 
on  slab  or  bar  samples.  The  usual  cutter  speed  used  in  propellant  milling  is 
approximately  900  rpm  and  the  table  feed  rates  are  varied  between  one  and 
ten  inches  per  minute,  depending  on  the  milling  characte ristic s  of  the  propel¬ 
lant.  The  general  philosophy  applied  is  that  the  best  specimens  are  produced 
by  a  combination  of  high  cutter  speeds  and  maximum  feed  rates. 

The  two  basic  cuts  milled  on  propellant  specimens  are  the  planar  surface  cut 
and  the  re-entrant  section  cut  on  a  JANAF  tensile  specimen.  A  simple,  manu¬ 
ally  operated  clamping  device  is  used  ordinarily  for  the  shaping  of  the  re¬ 
entrant  sectionof  the  JANAF  samples.  Baror  slab  samples  for  milling  planar 
surfaces  can  be  clamped  by  manual  vises  (Figure  A-5);  however,  this  tech¬ 
nique  has  restrictions  imposed  by  difficulties  in  maintaining  tolerances  during 
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Figure  A -3  Milling  Machine  Modified  for  Propellant  Machining 
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Figure  A-5  Propellant  Specimens  Being  Machined  in  Manual 
Vise  Mill  Table  Attachment 


milling,  particularly  when  slabs  or  long  bars  are  being  processed.  To  over¬ 
come  this  difficulty,  a  vacuum  gripping  plate  vise  may  be  utilized.  A  vacuum 
plate  vise  in  use  for  milling  tab  end  specimens  is  shown  in  Figure  A-6. 
A  standard  portable  vacuum  pump  maintains  satisfactory  gripping  of  speci¬ 
mens  . 


A- 1.4  PREPARATION  OF  JANAF  TENSILE  SPECIMENS 


JANAF  tensile  specimens  are  prepared  from  rough  cut  blanks  approximately 
five  inches  long  by  five -eighths -inch  thick  by  one  and  one -eighth -inch  wide. 
The  specimens  are  processed  in  groups  of  five  through  the  milling  operation, 
maintaining  in  the  final  specimen  the  tolerances  and  dimensions  shown  in 
Figure  A-7.  'fhe  surface  roughness  of  all  milled  specimens  is  considerably 
below  one  mil . 

In  the  preparation  of  JANAF  specimens,  tolerances  of  4  mils  are  maintained 
on  the  re-entrant  section  of  the  specimen. 


A -1.5  PREPARATION  OF  WOODEN  TAB  END  TENSILE  SPECIMENS 


Wooden  tab-end  specimens  are  prepared  using  case -bonding  techniques  in 
which  the  tab  segments  of  the  specimens  are  bonded  to  the  propellant  by  cast¬ 
ing  uncured  propellant  against  the  treated  wooden  surface.  The  wood  ends  of 
the  containers  are  cross  cut  into  one -half  by  six-by-six  inch  slabs  that  are 
used  to  make  box  ends  spaced  five  inches  apart.  Prior  to  casting  propellant, 
the  box  is  lined  with  a  suitable  bonding  agent.  In  practice  (PBAA,  Polysulfide) 
the  best  bonding  agent  has  proven  to  be  the  motor  case  liner  material.  After 
lining,  the  box  is  vacuum  cast  with  propellant  (Figure  A-8). 

Boxes  with  cured  propellant  are  sawed  into  slabs  of  the  desired  size,  retain¬ 
ing  the  wooden  sides  of  the  box  as  integral  parts  of  the  slabs.  The  sawed 
slabs  are  milled  on  both  surfaces  to  the  desired  thickness  and  sawed  into 
strips.  The  sawed  surface s  of  the  strips  finally  are  milled  to  obtain  rectangu¬ 
lar  bars  of  square  cross-section  (0,375  or  0.500-inch  square),  maintaining 
tolerances  of  2  mils,  The  specimen  is  illustrated  in  Figure  A -9- 

When  marginal  bonding  is  obtained,  wooden  tab-end  specimens  satisfactory 
for  use  in  viscoelastic  characterization  may  be  prepared  in  the  form  of 
JANAF/tab-end  hybrids  as  shown  in  Figure  A- 10.  These  specimens  display 
a  constant  gage  length  behavior  to  a  satisfactory  accuracy,  but  require  speci¬ 
alized  techniques  for  machining. 
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Figure  A-6  Tab  End  Tensile  Specimens  Being  Machined  on  Vacuum 
Vise  Mill  Table  Attachment 
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Figure  A -8 


Tab  End  Tensile  Sample  Box  Cast  with  Propellant 
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Figure  A- 10  JANAF  Tab  End  Hybrid  Tensile  Specimen 
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A- 1.6  PREPARATION  OF  SPECIMENS  FOR  BIAXIAL  TESTING 


Biaxial  propellant  specimens  are  prepared  by  casting  long  cylinders  (in  "Seal - 
rite"  containers,  for  example)  of  propellant  with  the  desired  diameter,  then 
cutting  and  milling  to  produce  coplanar  faces  (Figure  A-ll.)  Extremely 
toxic  or  expensive  materials  are  cast  directlyinto  molds  of  the  desired  shape 
and  tested  after  removal  from  the  mold  without  further  preparation. 


A-1.7  TENSILE  SPECIMEN  STRAIN  MEASUREMENT 


Two  methods  of  measuring  the  strain  in  JANAF  specimens  are  now  in  use  in 
the  solid  propellant  industry.  The  first  method  consists  of  placing  thin  lines 
or  bench  marks  on  the  central  portion  of  a  JANAF  specimen  and  photographi¬ 
cally  recording  their  separation  during  tests.  Data  reduction  for  this  process 
is  expensive  and  time  consuming.  The  method  in  use  at  LPC  (LPC-556-F-  1 , 
1963)  utilizes  a  plastic  extensometer  fixed  directly  on  the  test  specimen 
(Figure  1-8,  Section  1).  The  plastic  extensometer  consists  of  a  positive  film 
transparency  (Figure  A-12),  which  is  affixed  to  a  specimen  at  Point  A  and 
in  alignment  with  an  indicial  bench  mark  on  the  unstrained  specimen  (Point 
B).  A 8  the  specimen  is  stretched,  the  increase  in  length  can  be  measured  by 
noting  the  passage  of  the  specimen  bench  mark  past  the  indicial  marks  C  on 
the  extensometer.  When  the  spacing  of  the  marks  C  are  as  shown  (Figure 
A-12)and  the  length  between  the  pointof  attachment  of  the  gage  and  the  speci¬ 
men  bench  mark  has  increased  successively  to  1.025,  1.05,  1.075,  1.10,1.15, 
etc.,  the  corresponding  strains  will  be  2-1/2,  5,  7-1/2,  10,  15,  etc.,  per¬ 
cent.  The  strain  readings  are  placed  directly  on  the  test  record  by  use  of  the 
tester  "Pip  marker”,  reading  the  strains  visually  during  the  test.  A  typical 
test  record  is  shown  in  Figure  A-13.  Strain  values  are  then  scaled  directly 
from  the  test  record  during  data  reduction. 

In  practice,  the  gage  is  assembled  to  the  propellant  samples  by  cementing  a 
0.020-inch  diameter  metal  wire  to  the  plastic  extensometer  at  point  A,  which 
is  subsequently  cemented  to  the  propellant  with  a  quick  drying  adhesive.  To 
keep  the  free  end  of  the  extensometer  from  curling  away  from  the  propellant, 
a  staple  is  placed  loosely  over  the  lower  edge  of  plastic  in  the  flow  region  of 
the  propellant  (Figure  1-8,  Section  1). 
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Figure  A- 12  LPC  Plastic  Extensometer 
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Section  A-2 


TEST  SPECIMEN  HANDLING  AND  STORAGE 


Absorbed  moisture  is  probably  the  single  most  important  environmental  con¬ 
taminant  that  must  be  controlled  during  handling  and  storage.  Recommended 
standard  procedure  in  this  connection  is  to  store  all  samples  and  specimens 
in  plastic  bags  with  laboratory  desiccant  (e.g.  ,  silica  gel  or  "Drierite")  from 
the  completion  of  cure  until  a  specimen  is  fixed  in  the  tested  jaws.  Exposure 
of  propellant  during  propellant  sawing  and  machining  is  minimized  and  pre¬ 
ferably  carried  out  in  an  air  conditioned  and  humidity  controlled  environment 
(77°F,  50#R.H.  or  less). 

Du  ring  handling  and  storage,  test  specimens  are  packed  flat  on  planar  surfaces 
(e.g.,  stiff  cardboard)  to  prevent  specimen  distortion.  Any  specimen  sub¬ 
jected  to  accidental  flexing  should  be  discarded.  Temporary  sample  or  speci¬ 
men  storage  is  preferably  done  at  moderate  temperature  (70°F  ±  10°F),  but 
is  not  a  strict  necessity. 
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Appendix  B 

OPERATOR  EQUATION  -  PRONY  SERIES  COEFFICIENTS1 


Standard  notation  has  been  adopted  in  this  report  for  the  differential  operator 
equation  and  the  associated  Prony  Series  form.  This  notation  and  basic  re¬ 
lations  between  coefficients  of  the  equations  are  presented  herein. 


1.0  The  basic  operator  equation  is  given  by 
P<r  (t)  =  Qe(t) 

where 


P  =  £  a. 

i=o  1 


ji 

at1 


(B-l) 


(B-2) 


and 


Q  =  X  b 
i=o 


a1 

1  at1 


(B-3) 


2.0  Defining  Pand  Q  in  terms  of  the  associated  transform  polynomials  gives 
n  ,  n 

T*  =  £  a.  Sl  =  a  .n  (S  +  a.)  (B-4) 

l  n  t  ■■  i  t'  '  ' 


1  =  0 

n 

£ 

i=o 


<5  =  £  b.  Sl  =  b  .n,  (S  +  0.) 

.  _  i  n  1=1  '  i' 


(B-5) 


3.0  As  used  in  the  text  of  this  report,  the  partial  fraction  residuals  contained 
in  the  Prony  Series  coefficients  for  equations  of  the  form 


* <*>  ■  Co  +  C„o*  *  C, 

1=1 


or 


c(t)  =  +  D_  t  +.  £  D 


-v 


-v 


(B-6) 


oo 


i=i 


'i 


(B-7) 


1  Reference:  LPC  556-F-l. 
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are  given  by 
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(B-9) 


where  the  are  arbitrary  constants  defined  by  the  stress  or  strain  input 
conditions. 

It  should  be  noted  that  the  A  residuals  are  normalized  with  an  =  1  and  the  B 
residuals  with  bn  =  1.  The  coefficients  with  subscript  "oo"  are  associated 
with  constant  rate  stress  or  strain  inputs  only  (LPC  556-F-l). 

4.0  The  relationships  between  the  Prony  Series  coefficients  “i-  Ai>  or  ei- 
Bi  and  the  operator  coefficients  a^  bi,  are  as  follows.  For  strain  deter¬ 
mined  functions,  the  coefficients  are  obtained  directly  and  the  product 

n 

n  (s  +  «i)  (B- 10) 
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can  be  formed.  The  coefficients  aj  of  the  resulting  polynominal  are,  by  defi¬ 
nition,  the  coefficients  of  the  P  operator  when  an  =  1 .  For  stress  determined 
functions,  the  coefficients  8j  are  obtained  directly  and  the  product 


n 

n  (s  +  e.) 

ui 


(B-ll) 


yields  the  Q  operator  coefficients  with  bn  =  1 . 

In  the  usual  case,  it  is  necessary  to  obtain  both  the  aj  and  bj  operator  coef¬ 
ficients  from  either  a  strain  determined  or  stress  determined  Prony  Series. 
For  strain  letermined  functions ,  the  ajare  obtained  from  Equation  B  -  1 0  above. 

The  b.  are  obtained  as  follows. 
i 


b  =  A  +  Ai  +  Az  +  ■  •  •  +  A 
no  ‘  n 


bn_(  =  Aq  M"  +  A,Mr  (»i  =  0)  +  A2m”  (<y2=0)  +  >  •  •  +  A^M1}  (»n=0) 

I 

!  (B- 1 2) 

bl  =  Ao  MH-.  +  Al  Mn-i  (®,  =  0)  +  A*  MU-i  (®2=°)  +  •  •  '  +  An  Mn-,  (V0) 


b 

o 


A 

o 


M 


n 

n 


where  is  defined  as  the  summation  of  all  possible  products  of  the  n 
coefficients  taken  r  at  a  time.  For  example 


V.  -  Ac  (or,  +  <%  +  ...  +  orj  +  Aj  (w2  +  a,  +  ■  •  .  +  o»n)  + 

(B- 1  3) 

A,  («i  +  «,  +  •••+  or  )+•••+  A  (»,  +  o2  +  •  •  •  +  a  ) 


For  stress  determined  functions,  the  b.  are  found  from  Equation  B-ll  above. 

The  a.  are  obtained  as  follows  1 

i 


a  =  B  +  Bi  +  B2  ■+•••+  B 
no  n 


a«  =  M1}  (01  =  0)  +  B$  M*i  ( 0* = 0 )  +  •  •  •  +  B  M1}  (S  =0) 

n-i  o  n  n 

t 

(B-14) 

bi  =  Bo  +  B,  (0,  =  O)  +  B2  (MO)  +  •  •  •  +  Bn  (Bn=0) 


bo  -  Bo  m; 
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Appendix  C 

UNIAXIAL  TENSILE  CREEP  COMPLIANCE  TEST 


Scope:  A  method  is  described  for  laboratory  measurement  of  the  creep  com¬ 
pliance  D(t)  under  constant  stress  loading  using  time -temperature  reduced 
variable  techniques. 

Test  Equipment:  In  general,  standard  laboratory  equipment  is  adequate.  A 
precision  temperature  controlled  chamber  is  required  with  provision  for  con¬ 
trol  of  humidity  in  the  chamber.  A  typical  experimental  set-up  is  described 
in  connection  with  description  of  the  test  procedure. 

Test  Specimen:  Bonded  tab-end  tensile  specimens  as  described  inAppendixA. 

Test  Procedure:  Creep  compliance  D(t)  under  isothermal  conditions  is  deter - 
mined  by  observation  of  the  strain  versus  time  e(t)  under  constant  stress 
oQ  loading: 


D(t)  = 


Since  the  cross-sectional  area  of  a  tensile  specimen  decreases  on  extension, 
it  is  necessary  to  provide  a  means  for  continuous  decrease  of  the  applied  load 
F  on  the  specimen.  This  is  accomplished,  for  example,  by  a  compensating 
spring  in  series  with  the  creep  specimen  (Figure  C-l). 


With  reference  to  Figure  C-l,  a  typical  test  set-up  consists  of  a  specimen 
connected  in  series  through  a  dial  displacement  gage  (with  spring  action  re¬ 
moved)  with  a  weight  connected  with  an  extended  spring  (Figure  C-2).  The 
test  is  initiated  with  the  specimen  under  initial  load,  but  with  the  spring  ex¬ 
tended  to  give  an  initial  force  Fz  acting  on  the  weight.  The  spring  is  selected, 
in  consideration  with  the  weight  W  (W  includes  the  weight  of  the  connecting 
elements)  to  have  a  spring  constant  k  determined  by 


k  = 


(lb/in.) 


where 

AQ  =  unstrain  cross-sectional  area  of  the  test  specimen  (in.z) 
1Q  =  specimen  gage  length  (in.) 
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Figure  C -2 


Constant  Stress  Creep  Compliance  Apparatus 
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Upon  release  of  the  weight  to  act  on  the  specimen,  the  stress  oQ  acting  on  the 
specimen  will  be 

ao  =  "E”  +  k  Als 

where  A ls  is  the  initial  spring  extension.  Defining  strain  in  the  specimen  as 

e  (t)  =  (in./ in.) 

o 

where  A  l(t)  is  the  increase  in  specimen  length  with  time  as  measured  by  the 
dial  gage;  the  measured  creep  compliance  is  then  given  by 


D(t)  = 


e(t) 

°o 


The  teat  can  be  initiated  by  either  auxiliary  automatic  equipment  or  simply  by 
gently  releasing  the  weight  by  hand. 

Data  Reduction:  The  time  necessary  for  the  weight  to  be  fully  supported  by 
the  specimen,  t,,  is  noted  and  the  compliance  D(t)  under  constant  tempera¬ 
ture  conditions  is  calculated  for  times  in  excess  of  10  t,.  Data  for  various 
temperatures  may  be  obtained  and  reduced  to  master  curve  form  by  reduced 
variable  techniques  as  illustrated  in  Sections  1.  1  and  3.  1.  Typical  creep 
compliance  data  at  70°F  are  illustrated  in  Figure  C-3. 
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Figure  C-3  Typical  Creep  Compliance  Data 
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Figure  D-3  Nomenclature  for  Star  Valley  Stress  Concentration 
Relationships,  Figures  D-4  through  D-8 
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Three -Slotted  Grains 


Figure  D-5  Star  Valley  Stress  Concentration  Factor  H  for 
Four -Slotted  Grains 
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Figure  D-6  Star  Valley  Stress  Concentration  Factor  H  for 
Five -Slotted  Grains 
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Figure  D-7  Star  Valley  Stress  Concentration  Factor  H  for 
Six-Slotted  Grains 
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Figure  D-8  Star  Valley  Stress  Concentration  Factor  H  for 
Eight-Slotted  Grains 
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Figure  D-9  Nomenclature  for  Grain  Slump  Relationships, 
Figures  D-10  through  D-14 
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Figure  D-10  Dimensionless  Axial  Displacement  of  Free 
Corner  of  Free-Free  Grain 
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Figure  D-12  Dimensionless  Axial  Displacement  U  of  Free 
Corner  of  Fixed-Free  Grain  Z2 
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